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ABSTRACT 



Modern naval warfare has been increasingly dependent upon the 
acoustic silencing of the participants. Constrained viscoelastic layer 
damping of vibrating elements is one method which can be used to meet 
acoustic silencing goals. This paper considers constrained viscoelastic 
layer damping treatments applied to a thick aluminum plate, including 
single layer, double layer, a milled pocket plate, and a milled “floating 
element” configuration. High modal damping values were obtained for 
each damping configuration. The Modal Strain Energy method, using 
finite element analysis to estimate modal loss factors, was investigated for 
use as a tool in constrained viscoelastic layer damping design. A 
comparison of experimentally measured frequency response and modal 
loss factors with those predicted by the modal strain energy method is 
presented to confirm the possible use of the modal strain energy method as 
a design tool. 
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L EMTRODUCTION 



Modem naval warfare, especially undersea warfare, depends heavily 
on the vessel being acoustically silent. A major source of radiated noise is 
the vibration of shipboard components. The reduction of these vibrations is 
of utmost importance if a ship is to accomplish its mission. One method of 
vibration damping that shows promise in damping over a broad spectrum 
of low frequency vibration is constrained viscoelastic layer damping. The 
constrained viscoelastic layer method uses the high energy dissipation 
characteristic of viscoelastic materials during periodic motion of shear 
deformation to absorb and dissipate the vibrational energy of the system in 
question. Unfortunately the design and analysis of such constrained 
viscoelastic layer systems is difficult, due in part to the following: 

• The material properties of vdscoelastic damping materials vary 
greatly with temperature and frequency. 

• Closed form solutions to the equations of motion for constrained 
layer system exist only for beams and plates with simple 
boundary conditions. 

• The exact complex valued eigenvalue analysis for constrained 
viscoelastic layer damping systems using the finite element 
method requires large amounts of computer storage and CPU time. 
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Johnson and Kienholz developed the Modal Strain Energy (MSE) 
method which uses the ratio of strain energy for each mode shape to 
approximate the modal damping of a structure for a given constrained 
viscoelastic damping system [Ref. 1]. This method is very attractive due to 
its simple concept and very useful because it can be applied to any general 
cases with arbitrary shape by using the finite element method. However, 
its effectiveness compared with experiments were reported for only a few 
cases. Maurer examined the effectiveness of the MSE method for two 
damped plate configurations: 1) a simple sandwich configuration, and 2) 
a plate with a milled pocket with damping material inserted and a welded 
cover plate acting as a constraining layer in a previous work [Ref.2]. 
However, he could not verify for the milled pocket plate case since the cover 
plate warped and delaminated itself from the damping material during 
welding, resulting in negligible damping [Ref. 2]. 

For certain naval applications the components to be damped will be 
thick in construction and may be exposed to an unfriendly environment. 
Therefore, in this research the pocket plate configuration and a second 
milled plate using a “floating element” in conjunction with constrained 
layer damping are investigated. These damping treatments are compared 
with the simple sandwich type treatments consisting of single and double 
constraining layers. Simple plate geometries were used to facilitate the 
experimental and computational effort. Therefore, this paper addresses 
the experimental testing and analysis of four thick aluminum plates, each 
with a different constrained viscoelastic layer damping treatment. The 
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effectiveness of the MSE method is investigated for each different damping 
treatment by evaluating its accuracy in modal damping value prediction 
compared with experimental results and usefulness as a possible design 
tool. 
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II. THEORY 



A. VISCOELASTIC MATERIAL 

Viscoelastic materials of interest for general naval 
applications are polymeric compounds made up of long molecular chains. 
These molecular chains can be strongly, or weakly, linked together, 
depending on their chemical composition and processing. The damping 
characteristics of viscoelastic arise from the deformation and recovery of 
the polymer network. Material properties of a viscoelastic material vary 
with temperature and frequency. As such, the damping characteristics of 
a system will vary as its operating environment changes. [Ref. 3] 

Temperature will have the greatest effect on the material properties 
of damping materials [Ref. 3]. This effect is shown in Figure 2.1, where 
four distinct regions are observed. The lowest temperature region is the 
glassy region where the -material’s storage modulus is at its maximum 
value, and the loss factor is at a minimum. In the glassy region the 
modulus decreases slowly with temperature increase, whereas the loss 
factor increases rapidly with temperature. The second region is the 
transition region where the modulus decreases rapidly with increasing 
temperature and the loss factor reaches its peak value. The third region is 
the rubbery region where both the modulus and loss factor are at low values 
and show little variation with temperature. The fourth, and last, region is 
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the flow region and characterizes the behavior of some materials, mostly 
ceramics, at high temperatures. It should be noted that the transition 
region may vary in width from 20 °C up to a width of 200 °C. [Ref. 3] 

The effect of frequency on viscoelastic materials is not as great as that 
of temperature. The modulus of the viscoelastic always increases with 
increasing frequency. The loss factor will initially increase with frequency, 
then peak, and subsequently decrease as frequency increases. A plot of 
storage modulus and loss factor versus frequency is shown in Figure 2.2. It 
should be noted that this plot is over a range of approximately ten decades, 
and hence it becomes obvious that a temperature change of a couple degrees 
v/ill have a much greater effect on damping than a minor change in 
frequency. [Ref. 3] 

Linear viscoelastic materials behave in a hysteretic manner under 
cyclic excitation. Therefore, the mechanical behavior of a viscoelastic 
material during steady state vibration is best described by using a complex 
stiffness, k* [Ref. 3]. 

k* = k(l + ip) (2.1) 

where, 

r\ = material loss factor 

The use of a complex stiffness then leads to the use of a complex Young’s 
modulus and shear modulus [Ref. 3]. 
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E* = E(i + in) 
G* = G(i + in) 



( 2 . 2 ) 

(2.3) 



This concept of the complex modulus is used in subsequent analysis. 

Viscoelastic material properties are commonly displayed using a 
“reduced frequency nomogram.” The reduced frequency nomogram 
displays the variation of the viscoelastic material’s loss factor and modulus 
with temperature and frequency. The “reduced frequency”, fat. is an 

empirically determined function that accounts for the viscoelastic’s 
temperature and frequency dependence, and allows data for wide range of 
temperature and frequencies to be plotted on the same graph [Ref 4]. The 
reduced frequency nomogram for 3M ISD - 112 is shown in Figure 2.3. To 
find the loss factor and modulus using the nomogram, enter with the 
desired temperature and frequency. Follow the frequency line horizontally 
and the temperature line diagonally down the page until the two intersect. 
Then go vertically up or down to intersect the shear modulus or loss factor 
curves. Finally, read the value of the shear modulus or loss factor 
horizontally from the scale on the left [Ref. 4]. 



B. CONSTRAINED VISCOELASTIC LAYER DAMPING 

A simple constrained layer damping treatment consists of a base 
layer (the structure to be damped), a damping layer, and the constraining 
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layer. This coniigxiration is shown in Figure 2.4 with the thicknesses of the 
damping and constraining layers exaggerated for clarity. 

The physical mechanism of damping can be explained by referring to 
Figure 2.4. When the base layer is deformed in a mode of vibration, the 
surface away from the neutral axis elongates, stretching the viscoelastic 
material. The top layer, being a stiff elastic material, tends not to elongate, 
and thereby “constrains” the viscoelastic material. Consequently, the cyclic 
motions of vibration induce a cyclic shearing strain in the viscoelastic. 
This cyclic shearing strain, together with its associated hysteresis loop 
cause the vibrational energy to be dissipated as heat. For the constraining 
layer to be effective, its stiffness should not exceed that of the base layer. 
[Ref. 4 & 5] 



C. SYSTEM EQUATIONS OF MOTION 

Continuous systems, such as plates, possess distributed 
characteristics of mass, damping, and stiffness. Classical vibration 
analysis of such systems involves the formation of a mathematical model 
that discretizes the system into a finite number of components in order to 
approximate the total system. Such a formulation results in the following 
equation: 

[M] (x(t)) + [C] (x(t)j + [K] {x(t)} = (F(t)} (2.4) 
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where, 



[M] = system mass matrix 
[C] = system damping matrix 
[K] = system stiffness matrix 
(F(t)) = external excitation vector 
{x(t)) = displacement vector 

For an undamped system without excitation, the above equation reduces to 
the eigenvalue problem. 

[M]{x(t)l+[K]{x(t))=0 (2.5) 

This equation is then transformed to modal space using the linear 
transformation: 

{x(t)j =[ 0 ] (q(t)) (2.6) 



where, 

[0] = modal matrix 

{q(t)l = modal response vector 

Using this linear transformation the equation of motion can then be solved 
for the undamped modal frequencies and mode shapes. 

To solve for the frequency response of a damped system, the linear 
transformation is applied to equation (2.4): 

[M] [ 0 ] (q(t)) + [C] [ 0 ] (q(t)) + [K] [ 0 ] {q(t)} = {F(t)} (2.7) 
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Assuming that the damping matrix [C] is proportional to a linear 
combination of the stiffness matrix [K] and mass matrix [M], the damping 
matrix can then be diagonalized using the same linear transformation 
used to diagonalize [K] and [M] in equation (2.5) above. The diagonal terms 
of the damping matrix then become (TijCOj), where qj equals the modal loss 
factor and (Oj is the natural frequency of the i^^ mode [Ref. 1] . Using this 
approximate diagonal damping matrix results in a system of uncoupled 
modal equations of motion: 

qj(t) + qiCOi<i(t) + of qi(t) = f,(t) (2.8) 



where, 



qi(t) = modal acceleration of i^^ mode 
q(t) = modal velocity of i^^ mode 
qi(t) = modal displacement of i^^ mode 
qi = modal loss factor of i^^ mode 
coj = i^^ natural frequency 
fj(t) = modal force in i^^ mode 




Assuming that a sinusoidal excitation produces a sinusoidal response, 



(f(t)} = (f)ei‘^ {q(t)l = {Qlei^^ 



the response for the i^^ mode is then solved to be: 



(2.9) 
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( 2 . 10 ) 



cof - 0)2 + jCOiTliO) 

Subsequently, the response of the physical system can be found using: 

(x(t))=[({)]{Q}6)‘** (2.11) 

D. MODAL STRAIN ENERGY METHOD 

The equations of motion used to define the response of a system with 
viscoelastic materials need a complex eigenvalue analysis. However, the 
actual solution of these equations may be quite difficult. This is especially 
true when the system to be analyzed is comprised of materials whose 
properties vary with both temperature and frequency. Finite element 
techniques are generally used to compute the response of complicated 
systems. However, for the case of varying material properties many time 
consuming and costly runs must be made with the material properties 
changing at each frequency increment [Ref 1]. In addition to the costly 
analysis of a single design configuration, changes in design options, design 
requirements, of the search for an optimum design can make the expense 
of finite element analysis too great. The development of the Modal Strain 
Energy (MSE) method by Johnson and Kienholz, however, makes the finite 
element analysis of complex viscoelastically damped structures a viable 
option [Ref 1]. 



The modal strain energies can be obtained for finite element analysis, and 
are a standard output option of the NASTRAN finite element code [Ref. 7]. 

The modal frequency response of the structure is then calculated 
using the modal loss factors found in equation (2.12). When computing the 
modal frequency response of a damped structure, the modal properties in 
the system matrices are assumed to be constant. However, viscoelastic 
materials have storage moduli which are frequency dependent. To account 
for this frequency dependence, Johnson and Kienholz devised the following 
correction factor to be applied to the modal loss factors calculated in 
equation (2.12) [Ref 8]. 




(2.13) 



where. 



= corrected modal loss factor at the r^^ mode 



G 2 (fr) = viscoelastic shear modulus at the r*-^ modal 
frequency 



G2,ref = reference viscoelastic shear modulus used in the 
frequency response calculation 



storage module tC or loss factor t; 




Figure 2.1. Variation of Viscoelastic Material 

Properties with Temperature [Ref. 3]. 




Figure 2.2. Variation of Viscoelastic Material 
Properties with Frequency [Ref. 3]. 
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Figure 2.3. Temperature Frequency Nomogram for 3M ISD - 112. 
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Figure 2.4. Single Constrained Layer Configuration. 
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HL DESIGN OF DAMPED PLATES 



A . GENERAL SPECIMEN CONFIGURATIONS 

For experimental testing and analysis purposes, four different, yet 
related, constrained layer damping treatments were selected in addition to 
an undamped “reference” plate. Two of the damping configurations were 
simple sandwich treatments consisting of one and two viscoelastic layers 
respectively. Another treatment was the “pocket plate” which was 
previously investigated by Maurer [Ref. 2]. The pocket plate was made from 
a solid plate which was then milled to accept damping material and a cover 
plate. The final damping treatment was a “floating element” configuration 
consisting of a solid plate milled to accept a double layer of damping 
material and a welded cover plate. Section views of these damping 
configurations are shown in Figure 3.1. 

The purpose of the' pocket plate is to protect the viscoelastic material 
from materials such as oil and salt water, which may harm the 
viscoelastic. Since previous attempts at using a welded cover plate were 
unsuccessful due to the heat of welding causing a delamination between 
the damping material and cover plate [Ref. 2], it was decided to move the 
viscoelastic material from the welding point by recessing it into a shallow 
pocket of its own as shown in Figure 3.1. 



The “floating element” concept evolved because the welded cover plate 
of the pocket plate configuration does not produce the damping reaction that 
a true constraining layer would provide. If the cover plate is welded to the 
surrounding structure it cannot deform in bending as much as an 
unwelded constraining layer, thereby causing a reduction in the damping 
capability of the system. By using a piece of metal in the milled pocket with 
dimensions slightly smaller than the surrounding pocket, along with two 
layers of viscoelastic and a welded cover plate, a true constraining layer 
effect should be obtained. This configuration is shown in Figure 3.1. 

In order to approximate a possible system to be damped, a plate 
with large dimensions was selected. The dimensions of the plates used for 
the damping treatments are 114.3 cm (45 in) in length and 38.1 cm (15 in) 
in width. In addition, possible naval applications for this type of damping 
would probably consist of thick plate members. For this reason it was also 
decided that thick plates would be used for the damping treatments. All the 
base plates and constraining layers were made of a standard 6061— T6 
aluminum alloy. The design and selection of vdscoelastic layer, base layer, 
and constraining layer thicknesses is discussed in the following sections. 



R DESIGN OF THE SINGLE DAMPING LAYER CONFIGURATION 

In an attempt to approximate system loss factors and hence 
determine viscoelastic and constraining layer thicknesses for maximum 



damping, a method developed by Nashif [Ref. 9] based on an analysis of 
simple sandwich damping systems by Ross, Kerwin, and Ungar [Ref. 10] 
was used. The Ross - Kerwin - Ungar (RKU) equations are base of the 
analysis of the simple sandwich system shown in Figure 3.2. 

To find the loss factors of the damped system, the fiexural rigidity of 
the system must first be determined. For the above system, the flexural 
rigidity, El, is written as [Ref. 9]: 



El = Mil + Mil + Mii + EiHiD^ 
12 12 12 

+ E2H2(H2i - D)2+ EsHsiHai - Bf 



- D) + EsHaiHai - D) 



Hai - D 

1 + g 



(3.1) 



where. 



E2hJh2i -^] + g(E2H2H2i + E3H3H31) 

(3-2) 

EiHi + ^ + g(EiHi + E2H2 + E3H3) 



H3, = : Ha ^ JJ2 


(3.3) 




(3.4) 


g= 

E3H3H2K2 


(3.5) 


E = Young’s modulus 




G = Shear modulus 





I = second moment of area 
H = thickness of member 
K = wave number 



Subscripts refer to the layers labeled in Figure 3.1. No subscript refers to 
the composite system 



For a simply supported plate the wave numbers and modal 
frequencies are found using [Ref.9]: 



— Nfim'' 



EH^g, 

12(l-v2)Hp 



(3.6) 



where, 




a = semi-wave length of the plate 
b = semi— wave width of the plate 
V = Poisson’s ratio of the composite plate 
p = density of the composite plate 
gc = gravitational constant 



(3.7) 



To introduce damping into the equations it is necessary to use the 
complex modulus concept expressed in Section II. Substituting the 
appropriate expressions for the complex shear and Young’s modulus into 
equations (3.1), and assuming that damping in the base layer, (ri;|^), is 

small, and that the extensional stiffness of the damping layer is small 



(since E2 « E]^ and E2 « E3), the following expressions can be arrived at 
[Ref. 9 ]: 



ER3 = EiH? + EM + - P - 3)re 

c2 + 



( 3 . 8 ) 



EH^I = E3H3T13 + -J2— (a - P - 3)im 
c^ + d"^ 



( 3 . 9 ) 



where, 



a = gEiHiE 3 H 3 Hii{c(l-ri 2 ri 3 ) + d{^2 + Ha) + j[c(n 2 + hs)]} (3.10) 



p = E1H1E2H2H31 [c + dri2 + j(c7i2 - d)] 



(3.11) 



3 = 2gE3H3E3H3H3:H3i fl-2n2^3-nlM2,l3.n3-^lT,3) | ^ 3 ^ 2 , 

1+ j[c(2ri2 + Ti3 - h^a) - d(l - 2ri2Ti3 - h^jj] 



c = EiHi(l + g) + gE3H3(l - ri2Ti3) 
d = gEiHiTi 2 + gE 3 H 3 (ri 2 + TI3) 

j=^ 



(3.13) 

(3.14) 

(3.15) 



These equations were then applied to estimate the loss factors of the simple 
three— layer sandwich plate. The equations can be simplified by assuming 
that damping in the constraining layer (^3) is negligible [Ref. 9 ]. 



Since the boundary conditions for the plate used in this research 
(free-free-free-free) do not correspond to the simply supported conditions on 
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which equation (3.7) is based, modal frequencies for the free-free case were 
estimated using results from finite element analysis of a free-free plate and 
equation (3.6). The natural frequencies of an undamped plate with the 
dimensions previously given, and a thickness of 1.91 cm (0.75 in) were 
found using a normal mode extraction in NASTRAN. By substituting these 
modal frequencies into equation (3.6) and estimate of the wave parameter 
for each mode, K^, was obtained. Then, by using an iterative procedure 
outlined in Reference [9], modal loss factors were estimated for different 
layer thicknesses over a temperature range of 0.0 °C to 37.8 °C (30 °F to 
100 °F). 



The previous equations are easily programmed to compute loss 
factors for a wide variety of conditions. The variation of viscoelastic 
material properties with temperature and frequency was accounted for 
using a curve-fit to the reduced frequency nomogram developed by Drake 
[Ref. 11]. The material data for the following curve-fit equations is from the 
University of Dayton Research Institute [Ref. 12]. 



logio(M) = logio(ML) + 






(3.16) 



logio(ETA) = logio(ETAFROL) 

-I- A{SL + SH) -I- (SL - SH)( 1 - Vl A^)] 

log,o(FR) = log,„|F)-ij|Ii (3.18) 
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(3.19) 



, _logio{FR)-logio(FROL) 

C 

where, 

M = viscoelastic modulus 
ETA = viscoelastic material loss factor 
FR = reduced frequency (Hz) 

F = frequency (Hz) 

T = temperature (°F) 

and, 

T0 =40 °C(104 °F) 

FROM = 2.0x104 Hz 
MROM = 4.75x106 Pa (688.94 psi) 
n = 0.275 

ML = 6.0x104 Pa (8.7 psi) 

ETAFROL = 1.08 
SL = 0.45 
SH = -0.55 
FROL = 5000 Hz 
0 = 2.5 

In addition to the above constants, a Poisson’s ratio of 0.49 and a density of 
0.909 gram per cubic centimeter (0.035 Ibm/in^) was used for ISD-112 
[Ref.l2]. The following material properties were used for 6061-T6 
aluminum [Ref. 13]. 
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E = 70 GPadOxlO^ psi) 

V = 0.33 

p = 2.7 gm/cm3 (0.0968 Ibm/in^) 

Using the previous equations and material properties, a computer 
program was written to compute estimated modal loss factors and 
frequencies for a variety of base layer, viscoelastic layer, and constraining 
layer thicknesses. A listing of this program appears in Appendix A. 
Modal loss factors were computed for base layer thicknesses of 9.53 mm 
(0.375 in) to 19.05 mm (0.75 in) in 3.18 mm (0.125 in) increments. For each 
base layer thickness, the viscoelastic thickness was varied from 0.38 mm 
(0.015 in) to 1.52 mm (0.060 in) in 0.38 mm increments, and the 
constraining layer thickness was varied from 1.59 mm (0.0625 in) to 6.35 
mm (0.25 in) in 1.59 mm increments. In addition, loss factors were also 
computed for a viscoelastic thickness of 0.127 mm (0.005 in). From the 
results of the analysis, a carpet plot [Ref. 3], was made for each of the base 
layer conditions. The carpet plot reflects, for the the first mode, the 
maximum loss factor and its corresponding temperature for each 
viscoelastic layer/constraining layer thickness configuration. The carpet 
plot for a base layer thickness of 12.7 mm (0.50 in) is shown in Figure 3.3. 
Based on the carpet plots and a desire for maximum damping, as well as a 
system which could be moved easily, a base layer thickness of 12.7 mm(0.50 
in), a viscoelastic thickness of 0.38 mm (0.015 in), and a constraining layer 
thickness of 6.35 mm (0.25 in) was selected. The total system thickness was 
approximately 19.05 mm (0.75 in). This total system thickness would be 
maintained for all subsequent damping configurations. 



I 
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To maintain continuity between damping systems, the milled “pocket 
plate” was given a viscoelastic thickness of 0.38 mm and a cover plate 
thickness of 6.35 mm for a total system thickness of 19.05 mm (0.75 in). In 
order to keep the heat of welding away from the viscoelastic material, the 
ISD-112 was recessed into a shallow pocket as indicated previously in 
Figure 3.1, and as shown in the pocket plate system arrangement in 
Figure 3.4. Detail drawings of the pocket plate are shown in Appendix C. 



C. DESIGN OF THE DOUBLE DAMPING LAYER 

The design of the double layer damping system was accomplished 
using a modification of the RKU analysis used in the previous section. The 
RKU equations are used by working from the top layer of the damping 
system down towards the base layer. As shown in Figure 3.5, the H3' 
constraining layer along with the H2' viscoelastic layer are combined with 
the Hi’ layer to form a three— layer system. Using the RKU equations, the 
stiffness of this system is computed and considered to be the equivalent 
stiffness of the top three layers of the total constrained layer damping 
system. The top three layers were then considered as a single layer with 
the equivalent stiffness previously calculated, and the RKU equations were 
again applied to compute estimated modal loss factors for the entire double 
layer damping system [Ref.3]. 
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To maintain continuity among all the damping configurations, a 
base layer thickness of 12.77 mm (0.50 in) was chosen, and a total system 
thickness of 19.05 mm (0.75 in) was maintained. A design for high 
damping was then selected by computing modal loss factors for the 
constraining and damping layer thickness combinations shown in Table 
3.1. The estimated modal loss factors for the first mode of vibration in each 
configuration are plotted as shown in Figure 3.6. A listing of the program 
used to compute the loss factors is in Appendix B. 

From the data presented in Figure 3.6, viscoelastic thickness of 
0.38 mm (0.015 in) and constraining layer thicknesses of 3.18 mm (0.125 in) 
were selected for the double layer configuration. This particular 
configuration estimates high damping over a wider temperature range 
than the other thickness combinations. 



TABLE 3.1. THICKNESSES USED IN CALCULATION OF DOUBLE 
LAYER MODAL LOSS FACTORS. 



Confieuration 


1 


2 


4 


5 




H2 (mm) 


0.38 


0.38 


0.76 


0.38 


1.14 


Hi' (mm) 


3.18 


2.38 


2.38 


2.38 


2.38 


H2’ (mm) 


0.38 


0.38 


0.76 


1.14 


0.38 


H3’ (mm) 


3.18 


3.18 


2.38 


2.38 


2.38 



The milled “floating element” plate uses the same viscoelastic and 
constraining layer thicknesses as the simple two-layer configuration. In a 
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design similar to that of the pocket plate, the floating element and both 
layers of viscoelastic are recessed into a milled opening as shown 
previously in Figure 3.1 and further described in the floating element 
system configuration of Figure 3.7. 
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Single layer configuration 
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Double layer configuration 



Milled “pocket plate” 



Plate with “floating element” 



Figure 3.1. Four Damping Treatment Configurations. 



CONSTRAINING LAYER 



DAMPING LAYER 



BASE LAYER 




Figure 3.2. Elements of a Simple Sandwich Damping System. 
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Maximum System Loss Factor 
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Figure 3.3. Carpet Plot of Maximum Loss Factors for 
a Base Layer with HI = 12.7 mm. 




Figure 3.4. Arrangement of the Pocket Plate Configuration. 
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CONSTRAINING LAYER 
VISCOELASTIC 



CONSTRAINING LAYER 
VISCOELASTIC 



BASE STRUCTURE 



H3’ 

H2' 

HU 

H2 



HI 



Figure 3.5. Configuration of the Double Constrained 
Layer System. 
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Estimated Modal Loss Factors 



0.4 
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0.2 - 
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20 



Figure 3.6. 




120 



Modal Loss Factors for Double Constrained Layer 
Configurations. 
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Figure 3.7. Floating Element System Configuration. 
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IV. EXPERIMENTAL RESULTS 



A. TESTING ARRANGEMENT 

Experimental testing was performed on each of the four damping 
configurations and the undamped reference plate. In order to approximate 
the free-free-free-free boundary condition, each plate was suspended from 
the roof of the testing chamber using elastic cords as shown in Figure 4.1. 
All of the tests were performed in a temperature controlled environmental 
chamber which enabled temperatures to be maintained within ±1 °C. The 
primary component and user interface was the Hewlett-Packard (HP) 
3562A Dynamic Signal Analyzer (DSA). The HP-3562A was used to provide 
a swept sine signal to a vibration generator, and analyzed the returning 
data signals. The HP-3562A was used to compute the frequency response 
and coherence over a range of 50 Hz to 1050 Hz using the discrete Fourier 
Transform in the swept sine mode. Ten averages were performed at each 
data point using a frequency resolution of 625 mHz per step. The source 
level output to the vibration generator was set at 1.5 volts. 

A schematic of the testing apparatus is shown in Figure 4.2. Swept 
sine source signals were fed from the output jack of the HP-3562A to a 
Wilcoxon F3 vibration generator via the piezoelectric output of a Wilcoxon 
PA7C power amplifier. The vibration generator was mounted 73.48 cm 
(28.93 in) from one end, and 12.7 cm (5.0 in) from the front edge of each 
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specimen as shown in Figure 4.3. An integral force transducer was 
mounted in the base of the vibration generator to measure the force input to 
the plate. This force signal was then fed to input channel one of the DSA 
via a PCB 462-A charge amplifier. Plate accelerations were recorded at 
various points using a PCB 303A— 03 accelerometer as shown in Figure 4.3. 
Acceleration data was fed to input channel two of the DSA via a PCB 482A05 
power supply. Frequency response and coherence data was then recorded 
on disk for further analysis. 

Temperatures within the testing chamber were maintained using a 
NESLAB RTE— 8 refrigerated circulating bath which pumped fluid through 
a small heat exchanger in the testing chamber as shown in Figure 4.2. In 
order to accurately monitor the temperature of the plates, a small 
thermocouple was inserted in the base of each plate. 



B. TESTING PROCEDURE 

1 . Undamped reference plate 

An undamped, reference, frequency response measurement 
was made at a temperature of 15.6 °C (60 °F) to set a standard response by 
which to measure the effectiveness of the damping treatments. The 
undamped frequency response was recorded over a frequency remge of 
50 — 1050 Hz using a resolution of 625 mHz per point in the DSA. 
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2. Damped plate measurements 



Frequency response measurements of the damped plates were 
made at a temperature of 15.6 °C (60 °F) at several nodes on the plates in 
order to capture the damped response of as many modes as possible. A 
representation of these nodes is shown in Figure 4.3. Responses were 
recorded over a frequency range of 50 - 1050 Hz with a resolution of 
625 mHz per point in the DSA. Zoom measurements were also made to 
capture better data for certain modes. Modal loss factors were then 
estimated from the frequency response and coherence measurements 
using a curve-fitting technique described in Reference [14]. 



C. SINGLE DAMPING LAYER RESULTS 

The single damping layer treatment was tested at 4.44 °C (40 °F), 

15.6 °C (60 °F), and 26.7 °C (80 °F) so that the effects of temperature on the 
damping treatment could be determined. A plot of the single damping 
layer frequency response at 15.6 °C is shown in Figure 4.4. The single layer 
damping treatment resulted in high damping with modal loss factors 
ranging from 0.223 at 53.4 Hz to 0.091 at 876.6 Hz. Due to the coupling of 
modes, loss factors for all modes were not measured. The frequency 
response of the damped plate is characterized by a frequency shift to the left 
and a smoothing of the frequency response when compared to the 
undamped reference plate. The single layer treatment was especially 
effective at reducing the frequency response of a mode cluster between 650 
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and 950 Hz. The frequency band of this cluster was shifted approximately 
200 Hz with the amplitudes of the responses of the modes being 
dramatically reduced. The single layer treatment was also effective at 
reducing the amplitude of the response peaks over the entire spectrum of 
measurement. On average, the highest peaks of the frequency response in 
the undamped condition were reduced by 25 decibels, a reduction of 17.8 
times. 



The effect of temperature on the damping was quite pronounced as 
shown in Figure 4.5. As the testing temperature was decreased, the 
viscoelastic layer became stiffer and damping levels were increased. A 
comparative listing of the loss factors at different temperature is in Table 
4.1 and a plot of the modal loss factors is shown in Figure 4.6. Figure 4.5 
shows the trend of increased damping with temperature decrease, and a 
corresponding shift of modal frequencies to the right as the viscoelastic 
becomes stiffer. These changes are especially discernible at the lower 
frequencies. 
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TABLE 4.1. IVIEASURED MODAL LOSS FACTORS FOR THE 



SINGLE LAYER AT DIFFERENT TEMPERATURE. 



4.44 °C (40 °F) 15.6 °C (60 °F) 26.7 °C (80 °F) 



f(Hz) 


n 


f(Hz) 


64.3 


0.217 


53.7 


102.8 


0.223 


89.1 


157.4 


0.203 


138.9 


208.4 


0.183 


188.1 


327.5 


0.158 


301.6 


440.1 


0.144 


424.8 


474.9 


0.205 


476.8 


653.1 


0.079 


608.1 


768.0 


0.135 


644.2 


852.0 


0.093 


722.5 


919.0 


0.104 


817.6 


953.5 


0.073 


876.6 



21 


flHz} 


n 


0.223 


49.0 


0.117 


0.145 


83.5 


0.089 


0.172 


129.6 


0.082 


0.184 


177.9 


0.072 


0.120 


243.9 


0.062 


0.109 


409.4 


0.053 


0.111 


450.3 


0.066 


0.067 


552.9 


0.068 


0.069 


586.7 


0.052 


0.130 


629.3 


0.042 


0.081 


681.4 


0.096 


0.091 


790.1 


0.049 




858.2 


0.056 




978.1 


0.069 
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D. DOUBLE DAMPING LAYER RESULTS 



The double layer damping configuration was also tested at 4.44 °C, 
15.6 °C, and 26.7 °C. The frequency response of this configuration at 15.6 °C 
as compared to the undamped reference plate is shown in Figure 4.7. 
Damping in the double layer configuration is also high, with modal loss 
factors ranging form 0.301 at 53.3 Hz to 0.107 at 832.4 Hz. Due to modal 
coupling loss factors for all modes were not measured. The frequency 
response of the two-layer configuration is also characterized by a dramatic 
reduction in response amplitude and a frequency shift to the left. The peak 
undamped responses were reduced by an average of 27 decibels, or a 
reduction of 22.4 times from the reference condition. 

The effect of temperature on the double layer damping treatment is 
shown in Figure 4.8. As with the single layer case, damping in the double 
layer configuration increased with a decrease in temperature. This 
configuration also shows the shift of modal frequencies to the right as 
temperature decreases and the viscoelastic becomes stiffer. Modal loss 
factors for the double layer configuration are listed in Table 4.2, and are 
plotted for comparison in Figure 4.9. 

To compare the effectiveness of the single layer and double layer 
configurations their frequency responses at 15.6 °C are plotted in Figure 
4.10 with a plot comparing modal loss factors in Figure 4.11. The two 
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responses are quite similar although the double layer configuration does 
show an increase of approximately 22 percent in modal loss factor. 



TABLE 4J2. MEASURED MODAL LOSS FACTORS FOR THE 

DOUBLE LAYER AT DIFFERENT TEMPERATURE 



4.44 °C (40 °F) 15.6 °C (60 °F) 26.7 °C (80 °F) 



ffHz) 


n 


f(Hz) 


H 


f(Hz) 


n 


51.9 


0.273 


55.3 


0.301 


49.7 


0.202 


122.5 


0.224 


87.1 


0.215 


81.9 


0.188 


212.7 


0.187 


142.1 


0.217 


127.6 


0.156 


290.5 


0.197 


190.5 


0.212 


172.4 


0.117 


382.7 


0.174 


297.8 


0.139 


278.2 


0.097 


494.2 


0.144 


366.8 


0.154 


317.7 


0.077 


558.2 


0.169 


419.3 


0.125 


396.1 


0.070 


602.2 


0.198 


441.2 


0.100 


428.5 


0.077 


&42.4 


0.157 


618.7 


0.098 


602.0 


0.067 


752.8 


0.168 


680.9 


0.096 


638.4 


0.072 


815.4 


0.159 


715.5 


0.060 


738.0 


0.071 


896.3 


0.156 


832.4 


0.107 


802.6 


0.077 


968.3 


0.130 






846.1 


0.050 
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E. POOiET PLATE RESULTS 



The milled pocket plate was constructed as previously shown in 
Figure 3.3. The constraining layer, or cover plate, was welded in place 
using tack welds in an attempt to keep the damping material away from the 
heat of welding, and the cover plate from warping, instead of using a 
continuous weld bead as was previously attempted [Ref. 2]. The cover plate 
was welded to the base at the corners, at the midpoint of the short side and 
at three equally spaced locations along the long dimension as shown in 
Figure 4.12. Following welding the plate was tested to ensure that the 
viscoelastic had not been damaged by the heat of welding. 

The pocket plate was tested at 15.6 °C (60 °F) and the frequency 
response is shown in Figure 4.13. The response indicates that the 
viscoelastic layer was not damaged by welding and that good damping was 
attained. Modal loss factors ranged from 0.067 at 62.1 Hz to 0.090 at 923 Hz. 
Although damping is good, it is approximately half that of the single layer 
configuration. One reason for this is that the viscoelastic material does not 
completely cover the base structure. Another reason is the presence of the 
welded cover plate. Due to the welded conditions the cover plate cannot 
induce shear deformation in the viscoelastic layer as well as a true 
constraining layer, and therefore produces less damping than the single 
layer configuration. The effects of the welded cover plate are especially felt 
in modes below 300 Hz where the frequency response is quite peaked. The 
response curve becomes more rounded and the effects of the damping layer 
are seen as frequency increases. 
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The modal loss factors for the pocket plate are listed in Table 4.3 and 
are plotted in Figure 4.14. Even though the damping is less than the single 
layer, the plate is still adequately damped as shown in the frequency 
response plot in Figure 4.13. In this configuration the modal loss factors 
remained relatively constant throughout the testing spectrum. The 
increase in modal loss factor values above 800 Hz is due primarily to modal 
coupling, and the measured modal loss factors in this range are not 
reliable. 



TABLE 4.3 MEASURED MODAL LOSS FACTORS FOR THE 



POCKET PLATE AT 15.6^C. 



UEzi n 



f(Hz) n 



62.1 0.067 

93.4 0.042 

142.0 0.056 

194.9 0.060 

308.9 0.051 

437.5 0.041 

495.6 0.036 



565.3 0.043 

625.7 0.047 

643.2 0.049 

687.0 0.044 

841.9 0.056 

891.5 0.077 

923.0 0.090 
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F. 



FLOATING ELEMENT RESULTS 



The milled “floating element” configuration was constructed as 
previously shown in Figure 3.7. The cover plate was welded in a fashion 
similar to the pocket plate as shown in Figure 4.12. The center 
constraining layer, or floating element, was made slightly smaller than the 
surrounding structure thus allowing the floating element to act as a “true” 
constraining layer. 

The frequency response of the floating element configuration at 
15.6 °C (60 °F) is shown in Figure 4.15. The floating element is quite 
effective as the response shows a good reduction in peak modal response. 
Measured modal loss factors range from 0.089 at 66 Hz to 0.064 at 935 Hz. A 
listing of measured modal loss factors is in Table 4.4 and are plotted in 
Figure 4.16. As with the previous cases, the frequency response of the 
floating element configuration is characterized by a frequency shift to the 
left and a smoothing of the response as frequency increases. 

In a comparison of the pocket plate and floating element 
configurations, the frequency responses are plotted in Figure 4.17. A 
comparison of modal loss factors for the two configurations is shown in 
Figure 4.18. The two frequency response plots are similar, however, the 
frequency response of the floating element configuration is more rounded 
than that of the pocket plate. The major difference between the two 
configurations is seen in Figure 4.18. Modal loss factors for the floating 
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element show an average increase of 25 percent over those of the pocket 
plate. Reasons for this increase are the added constraining effect of the 
floating element and additional layer of damping material present. 



TABLE 4.4 MEASURED MODAL LOSS FACTORS FOR THE 
FLOATE'JG ELEMENT AT 15.6°C. 



f(Hz) 


n 




n 


66.0 


0.089 


508.1 


0.089 


104.9 


0.040 


540.6 


0.115 


154.0 


0.094 


618.7 


0.121 


211.4 


0.058 


656.7 


0.100 


272.8 


0.063 


708.6 


0.088 


323.5 


0.075 


857.6 


0.089 


435.8 


0.090 


935.4 


0.064 
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SUPPORTING BEAM 




Figure 4.1. Testing Configuration in Testing Chamber. 
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Figure 4.2. Schematic Diagram of Testing System. 
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Figure 4.3, Shaker and Accelerometer Locations. 
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Fxd Y 50 Hz SINGLE 15. S C 1.05k 



Figure 4.4. Frequency Response of the Single Layer 
Configuration at 15.6 °C. 



[ : undamped response; 



damped response] 




Fxd Y 50 Hz 4.44 / 15.6 / 26.7 C 1.05k 



Figure 4.5. Frequency Response of the Single Layer 

Configuration at Different Temperatures. 

[ : 4.44 °C, : 15.6 °C, : 26.7 °C] 
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Modal Loss Factors 




Fi^re 4.6. Comparison of Modal Loss Factors for the Single 
Layer Configuration at Different Temperatures. 




FxdXY 50 Hz DOUBLE LAYER 15.6 C 1 . 05k 



Figure 4.7. Frequency Response of the Double Layer at 15.6 °C. 
[ : undamped response , : damped response] 
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FxdXY 50 Hz 4 . 44 / 15 . 6 / 26 . 7 C 1 . 05k 



Figure 4.8. Frequency Response Comparison of the Double 
Layer Configuration at Different Temperatures. 

[ : 4.44 °C, : 15.6 °C , ; 26.7 °C] 
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Modal Loss Factors 
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Figure 4.9. Comparison of Modal Loss Factors for the Double 
Layer at DifTerent Temperatures. 
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FxdXY 50 Hz SINGLE / DOUBLE 1.05k 



Figure 4.10. Frequency Response of the Double and 
Single Layer at 15.6 °C. 
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Figure 4.11. Modal Loss Factors for Single and double layer 
Configurations. 
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Figure 4.12. Location of Tack Welds on the Cover Plate 
of the Pocket Plate Configuration. 
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Figure 4.13. Frequency Response of the Pocket Plate 
Configuration at 15.6 °C. 

[ : reference plate , : pocket plate ] 
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Figure 4.14. Modal Loss Factors for the Pocket Plate 
Configuration. 
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Figure 4.15. Frequency Response of the Floating Element 
at 15.6 °C . 

[ : undamped plate , : floating element ] 
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Figure 4.16. Modal Loss Factors for the Floating Element 
Configuration. 
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Figure 4.17. Comparison of Frequency Responses for the 

Pocket Plate and Floating Element Configurations. 
[ ; pocket plate , : floating element ] 
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V. FINITE ELEMENT RESULTS 



A . UNDAMPED REFERENCE PLATE 

The first step in the finite element analysis procedure was to model 
and analyze the undamped reference plate for its modal frequencies and 
frequency response. The finite element model was generated using 
PATRAN, a computer aided interactive graphics program developed by 
PDA Engineering. PATRAN is widely used for conceptual, preliminary, 
and detailed design and analysis of complex systems. One of PATRAN’s 
major advantages is in the interactive construction of finite element models 
for use by MSC/NASTRAN, and its ability to display MSC/NASTRAN 
results in an easily understood graphic format [Ref. 15]. 

The reference plate was modeled using 84 plate (QUAD4) elements as 
shown in Figure 5.1. The QUAD4 element is an isoparametric element 
with four nodes, one at each corner of the element [Ref. 7]. A normal mode 
extraction was then performed in order to compare numerical results with 
experimentally determined modal frequencies. Once satisfactory 
agreement between the modal frequencies calculated in NASTRAN and 
those obtained experimentally was attained, a direct frequency response 
calculation was performed in NASTRAN. This frequency response was 
then used as the reference response for further finite element models 
incorporating viscoelastic damping treatments. The excitation for the 
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frequency response calculation was a sinusoidal force with an amplitude of 
1.0 applied at node 66. The response point was node 58 as shown in Figure 
5.1. These two nodes correspond to the points on the plates used in the 
experimental portion of the research where the vibration generator and 
accelerometer were attached. 



B. SINGLE DAMPING LAYER 

The modeling of the single constrained layer damping system was 
done using techniques described by Johnson and Kienholz [Ref. 1]. As 
shown in Figure 5.2, the viscoelastic layer is modeled using solid (HEXA) 
elements, while the base layer and constraining layer were modeled using 
QUAD4 elements. The HEXA element is a solid, isoparametric element 
having eight nodes, one at each corner of the element with three 
translational degrees of freedom at each node [Ref. 7]. The use of solid 
elements for the viscoelastic layer allows the strain energy due to shearing 
to be adequately represented. Plate elements are used in the base layer and 
constraining layer because of their ability to account for stretching and 
bending deformations. The plate element allows its nodes to be offset from 
the plate’s center to the surface of the plate, coincident with the corner 
nodes of the solid viscoelastic elements [Ref. 1]. Thus, the single 
constrained layer system was modeled using only two layers of nodes, a 
simple process in PATRAN. For the single layer configuration, a model 
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having 84 elements per layer, and an element meshing scheme shown in 
Figure 5.2 was used. 

Once the damped plate had been modeled, normal mode extractions 
were made using MSC/NASTRAN. Five separate runs were conducted 
using the material properties of ISD— 112 at 50, 200, 500, 800, and 1000 Hz 
and at a temperature of 15.6 °C (60 °F). In addition to the modal 
frequencies, the strain energy in the viscoelastic elements and the entire 
model were output from NASTRAN. 

Since the shear modulus of a viscoelastic material changes with 
frequency, it was necessary to estimate the actual modal frequencies of the 
damped plate using an interpolation procedure outlined by Johnson and 
Kienholz [Ref. 16]. The first step of the interpolation process was to plot the 
shear modulus of ISD-112 versus frequency from 5 to 1000 Hz. Then, for 
the first mode, using NASTRAN results based on ISD-112 material 
properties a 50 Hz, the first modal frequency predicted by NASTRAN and 
the corresponding shear modulus were plotted. The same was then done 
using the first natural frequency predicted by normal mode extraction 
based on ISD-112 material properties at frequencies of 500 and 1000 Hz. A 
curve was then passed through these three points. The point where the 
NASTRAN modal frequencies for the first mode intersected the ISD-112 
shear modulus curve was taken to be the interpolated modal frequency of 
the single damping layer configuration. A plot of the intersection of these 
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two curves for the first and second modes is shown in Figure 5.3. This 
interpolation process was then repeated for each mode through 1000 Hz. 

Once the interpolated modal frequencies were found, the modal 
strain energy equations (2.12) and (2.13) were used to compute modal loss 
factors for the single layer configuration. A set of modal loss factors was 
computed based on the modal strain energies computed using viscoelastic 
properties at reference frequencies of 50 ,200, 500, and 800 Hz. A set of 
composite modal loss factors for the modal frequencies near these reference 
frequencies was then selected. The resulting modal loss factors are shown 
in Table 5.1 and are plotted versus frequency in Figure 5.4. As seen in 
Figure 5.4, the modal strain energy method is predicting high damping for 
this configuration with an average modal loss factor of 0.195. 



TABLE 5.1. ESTEVIATED MODAL LOSS FACTORS FOR THE 
SINGLE LAYER USING THE MODAL STRAIN 
ENERGY METHOD 



f(Hz) 


n 


f(Hz) 


n 


63 


0.249 


517 


0.237 


95 


0.181 


586 


0.217 


153 


0.234 


632 


0.164 


202 


0.212 


663 


0.155 


285 


0.229 


712 


0.180 


322 


0.214 


827 


0.134 


447 


0.200 


869 


0.152 


463 


0.108 


883 


0.115 


483 


0.258 
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Using the set of composite modal loss factors, the modal frequency 
response of the damped plate was computed using MSC/NASTRAN. Modal 
damping was introduced to the model using the SDAMP option in the Case 
Control Deck and the TABDMPl damping table in the Bulk Data Deck as 
described in the MSC/NASTRAN Handbook for Dynamic Analysis [Ref. 17]. 
Since NASTRAN uses a linear interpolation between points in the damping 
table to describe the modal damping in the model [Ref 7], a simple curve fit 
was applied to the set of composite modal loss factors as shown in 
Figure 5.4. Points from this curve fit were then used in the NASTRAN 
damping table. To compute the modal frequency response, a unit excitation 
force was applied at the same node as the undamped plate, and the node 
used for the response was also the same as the undamped plate. 

The results of the modal frequency response calculations are shown 
in Figure 5.5. The dashed line represents the undamped reference plate, 
and the solid line represents the modal frequency response of the single 
layer configuration. Material properties at 200 Hz were used for the 
ISD-112 damping material. The first thirty modes were used in the modal 
summation for the response. A listing of the MSC/NASTRAN data deck 
used to compute the modal frequency response is in Appendix D. 

The modal loss factors estimated using the modal strain energy 
method are compared to those measured experimentally for the single layer 
configuration in Figure 5.6. The estimated loss factors are greater than the 
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experimentally determined loss factors throughout the spectrum of 
interest, and especially in the middle frequencies. 

The modal frequency response of the single layer configuration is 
compared to the experimentally measured frequency response in 
Figure 5.7. The comparison was accomplished by normalizing both the 
experimentally determined frequency response and the frequency response 

computed in NASTRAN. Both responses were normalized using a value of 
in/s0C^ 

1.0 — — . The effects of the greater loss factors estimated by the modal 

strain energy method are obvious as the level of the predicted response is 
lower than the measured response. The shift in frequency between the two 
curves is due to the finite element model being inherently stiffer than the 
actual system. The correlation between the two curves is especially good 
below 250 Hz as this is where the differences between estimated and 
measured modal loss factors are the least. 



C. DOUBLE DAMPING LAYER 

The modeling of the double constrained layer damping system was 
accomplished as shown in Figure 5.8. The double layer configuration 
consists of a base layer modeled with offset QUAD4 elements, two 
viscoelastic layers consisting of HEXA elements, and the top constraining 
layer modeled with offset QUAD4 elements. The middle constraining layer 
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was modeled using three layers of HEXA elements in order to give this 
layer the stiffness necessary to act as a constraining layer. The model was 
meshed using 60 elements in each layer as shown in Figure 5.8. 

Using this model, the modal strain energy method was applied to 
determine approximate modal frequencies and loss factors for the double 
layer configuration. To determine the loss factors, normal mode 
extractions were performed using reference frequencies of 50, 200, 500, 800, 
and 1000 Hz. A composite set of modal loss factors for the double layer 
system was then compiled based on the estimated modal frequency’s 
relation to the reference frequency used to calculate modal strain energies. 
This composite set of modal loss factors is listed in Table 5.2, and is plotted 
in Figure 5.19 as a comparison to the experimentally determined modal 
loss factors for the double layer configuration. The estimated modal loss 
factors for the double layer show high damping, but they compare favorably 
with those measured experimentally. 

The modal frequency response of the double layer configuration was 
computed in a manner similar to the single layer in that smoothed loss 
factor data was used in the MSC/NASTRAN damping table. Likewise, a 
unit excitation force was applied, and the first thirty modes were used in 
the modal summation. The results of the modal frequency response 
calculation are shown in Figure 5.10. The dashed line represents the 
undamped response, and the solid line represents the frequency response of 
the double layer configuration. 
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TABLE 5.2. ESTIMATED MODAL LOSS FACTORS FOR THE 



DOUBLE LAYER USING THE MODAL STRAIN 
ENERGY METHOD. 



f(Hz) 


n 


f(Hz) 


n 


59 


0.278 


544 


0.209 


89.5 


0.219 


578 


0.141 


140 


0.239 


593 


0.131 


189 


0.213 


659 


0.177 


258 


0.214 


754 


0.102 


300 


0.205 


783 


0.076 


402 


0.179 


802 


0.144 


427 


0.176 


956 


0.077 


456 


0.264 


969 


0.111 


486 


0.242 


994 


0.048 



The frequency response calculated using NASTRAN was compared 
to the experimentally determined frequency response of the double layer 
configuration as shown in Figure 5.11. Comparison of the two frequency 
response curves shows similarity in form, but a much lower response level 
for the numerically determined response. Once again this could be due to 
the higher damping predicted by the modal strain energy method and the 
inherently higher stiffness of the finite element model. 
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D. POCKET PLATE RESULTS 



The pocket plate configuration was modeled with the same offset 
plate elements and solid viscoelastic elements as the single layer 
configuration. However, the pocket plate required the modeling of the 
milled structure around the viscoelastic material and the welds between 
the cover plate and milled plate. A representation of the model is shown in 
Figure 5.12. The base structure, cover plate, and the structure immediately 
around the cover plate was modeled using offset plate elements. The 
viscoelastic material and the portion of the structure immediately adjacent 
to it were modeled using the solid HEXA elements. Since the viscoelastic 
material and cover plate are physically separated from the surrounding 
plate, except where the viscoelastic is adhered to the base structure, care 
was necessary in creating the finite element mesh. 

The model was created in PATRAN using PATRAN’s node editing 
and equivalencing capabilities [Ref. 18]. This allowed the generation of a 
finite element mesh with two nodes at the same geometric point in space. 
Using this node editing technique, a mesh was created which allowed the 
viscoelastic and cover plate to vibrate separately from the surrounding 
structure, yet at the same time, keep the number of elements and nodes in 
the model to a minimum. The welded points on the cover plate were also 
modeled using node editing techniques. At weld points the finite element 
node on the cover plate was equivalenced with its corresponding node on the 
base structure, resulting in a single node and a connection between an 
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otherwise separate base structure and cover plate. At non— welded points 
on the cover plate there were two nodes at the same geometric point; one to 
represent the cover plate, and the other to represent the base structure. The 
model was meshed using a 5x11 mesh resulting in 40 elements in the cover 
plate, viscoelastic and base layer as shown in Figure 5.12. 

Using this modeling scheme, the modal strain energy method was 
employed to estimate the modal frequencies and loss factors. Normal mode 
extractions were made using viscoelastic material properties at 50, 200, 500, 
800, and 1000 Hz. Using these reference frequencies a composite set of 
modal loss factors was obtained. These loss factors are listed in Table 5.3 
and are plotted versus frequency in Figure 5.13. The estimated loss factors 
give good damping over the spectrum of interest with an average modal loss 
factor of 0.075. As with the pre\dous cases, damping values used for the 
modal frequency response calculation came from a curve— fit to the set of 
composite modal loss factors. 

The modal frequency response of the pocket plate was computed 
using the first 30 modes and viscoelastic material properties at 200 Hz. The 
resulting estimated frequency response is shown in Figure 5.14. The 
response shows a definite frequency shift to the left along with good 
damping of the frequency response when compared to the undamped 
response. 
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The estimated modal loss factors and modal frequency response for 
the pocket plate were compared to those measured experimentally. The 
loss factor comparison is shown is Figure 5.15 and the frequency response 
comparison is shown in Figure 5.16. The estimated modal loss factors are 
higher than those measured experimentally, however, the frequency 
responses compare quite favorably with each other. The frequency 
response curve calculated through finite element analysis has a lower 
response level and a frequency shift to the right of the measured frequency 
response. This is expected due to the increase in damping predicted by the 
modal strain energ)^ method and by the fact that the finite element model is 
inherently stiffer than the physical system 



TABLE 5.3. ESTIMATED MODAL LOSS FACTORS FOR THE 



POCKET PLATE USING THE MODAL STRAIN 



ENERGY METHOD. 



f(Hz) Q 



flH^ n 



67 0.046 

109 0.034 

162 0.091 

217.5 0.063 

295 0.091 

333 0.047 

453 0.07 

477 0.093 



511 0.113 

588 0.104 

625 0.086 

668 0.089 

713 0.077 

831 0.077 

834 0.065 

868 0.051 
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Figure 5.1. Finite Element Model of the Undamped 
Reference Plate. 
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DIMENSIONS IN CENTIMETERS (INCHES) 



Figure 5.2. Finite Element Representation of the 

Single Constrained Layer Configuration. 
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Figure 5.3. Intei-polation of the First and Second Modal Frequencies 
for the Single Constrained Layer Configuration. 
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Modal Loss Factor 




Figure 5.4. Estimated Modal Loss Factors for the Single Layer 
Configuration with the Curve Fit Used for the 
MSC/NASTRAN Damping Table. 
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Figure 5.5. Calculated Modal Frequency Response of the Single 
Layer Configuration Using NASTRAN. 

[ : reference plate , : single layer] 
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Figure 5.6. Comparison of Estimated and Measured Modal Loss 
Factors for the Single Layer Configuration. 
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Figure 5.7. Comparison of Estimated and Measured Frequency 
Response for the Single Layer Configuration. 
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Figure 5.8. Finite Element Representation of the Double Layer 
Configuration. 
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Figure 5.9. Modal Loss Factors for the Double Layer Configuration 
as Determined from the Modal Strain Energy Method 
and Determined Experimentally. 
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Figure 5.10. Calculated Modal Frequency Response for the Double 
Layer Configuration Using NASTRAN. 

[ ; reference plate , : double layer] 
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Figure 5.11. Comparison of the Experimentally Determined 
and Numerically Predicted Frequency Responses 
for the Double Layer Configuration. 
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Figure 5.12. Finite Element Representation of the 
Pocket Plate Configuration. 



84 



Modal Loss Factor 



0.12 



0.10 
0.08 
0.06 
0.04 
0.02 

0 200 400 600 800 1000 

Frequency (Hz) 

Figure 5.13. Estimated Modal Loss Factors 
for the Pocket Plate. 
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Figure 5.14. Estimated Frequency Response 

for the Pocket Plate Configuration. 

[ : reference plate , : pocket plate] 
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Figure 5.15 Comparison of Experimentally Measured and 
Estimated Modal Loss Factors for the Pocket 
Plate Configuration. 
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Figure 5.16. Comparison of Experimental and Predicted 
Frequency Responses for the Pocket Plate 
Configuration. 
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VI. CONCLUSIONS 



Constrained viscoelastic layer damping is an extremely effective 
method for reducing broadband vibration. In each of the experimental 
cases the peak amplitudes of frequency response were reduced by 
approximately 25 decibels, a reduction of 18 times below the undamped 
reference plate. 

In a comparison of experimentally determined modal loss factors for 
the four treatments, the double layer configuration yields the largest 
damping and the pocket plate the least as shown in Figure 6.1. Of 
particular note are the performance of the pocket plate and floating element 
configurations. Although they are not ideal configurations in terms of 
“true” constrained viscoelastic layer damping, the damping levels achieved 
are quite satisfactory. As was previously reported in Section IV, the 
floating element configuration yielded an average increase of 25 percent 
over the modal loss factors of the pocket plate. It is also noted that the 
average modal loss factor for the pocket plate is approximately 50 percent of 
the average modal loss factor of the single layer configuration. Similarly, 
the average modal loss factor for the floating element configuration is 
approximately 50 percent that of the double layer treatment. 

The modal strain energy method tends to overpredict the modal loss 
factors by as much as 50 percent for these highly damped, thick plates. 
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Figure 6.2 shows that, although the modal strain energy method predicted 
modal loss factors greater than those measured, the same relative 
differences in modal loss factor between damping configurations are 
maintained. The estimated modal loss factors for the pocket plate are 
approximately 60 percent less than those of the single layer configuration. 
This indicates that there is a consistency between the damping values 
predicted by the modal strain energy method and those of the actual 
physical system. 

There are several possible reasons for the differences between the 
experimentally determined modal loss factors and those estimated by the 
modal strain energy method. The first is that the material properties of 
ISD-112 as reported by the 3M Corporation on the reduced frequency 
nomogram may not be consistent with the material properties actually 
present in the material used. Since the numerical analysis was based on 
the reported material properties this is a possible source of uncertainty in 
the results. 

The second source of uncertainty is in the adhesion of the ISD— 112 to 
the aluminum plates. Although the plates were clean when the ISD— 112 
was applied, it was noted that the adhesive qualities of the ISD— 112 were not 
uniform throughout the material. A lack of adhesion may cause a 
reduction in the damping capability of the system. 
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Another source for the difference between the experimental results 
and numerical results lies within the finite element model. The number of 
elements used in the model greatly affects its “stiffness.” As the number of 
elements is increased the model should become less stiff and results are 
expected to approach those that are measured. Also, the type of element 
used to model the base layer and constraining layer may have an effect. In 
this research plate elements were used to model both the base and 
constraining layers. It is possible that one or more layers of solid elements 
may produce results that agree better with experimental results. 

One drawback to the modal strain energy method in design is the 
large amount of CPU time required for normal mode extraction and modal 
frequency response; especially in complex structures with a large number 
of elements. Therefore, although the modal strain energy method is good 
for analyzing a design, it may face a big difficulty to be used for design 
optimization due to the large amount of CPU time required for normal 
mode extraction and modal frequency response calculations. 
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Figure 6.1. Comparison of Experimentally Measured Modal Loss 
Factors for the Four Damping Configurations. 
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Modal Loss Factor 




Figure 6.2. Comparison of Estimated Modal Loss Factors for 
the Single Layer, Double Layer, and Pocket Plate 
Configurations. 
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VII. RECOMMENDATIONS 



The thick plate used in this research is a generic model of many 
physical systems that may see use in naval application. There are several 
areas which deserve more research and clarification, including the 
following: 

• Model the floating element configuration in finite elements and 
check the effectiveness of the modal strain energy method in 
predicting modal loss factors for this configuration. 

• Investigate the relationship between mode shapes and damping 
values. It seems some modal damping values are very low due to 
their twisting mode shapes. 

• For pocket plate and floating element configurations, weld the whole 
cover plate (continuous weld) and investigate the damping 
characteristics to compare the vibration reduction capabilities of 
the two different milled plate damping treatments. 

• Investigate methods for improving the adhesion of ISD-112 to the 
base structure and constraining layer. 
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APPENDIX A 



FORTRAN PROGRAM USED TO COMPUTE MODAL LOSS FACTORS 
FOR THE SINGLE DAMPING LAYER DESIGN 

This program used the Ross-Kerwin-Ungar equations of Section III 
to compute estimated modal frequencies and loss factors for the single layer 
configuration. Modal frequencies for an undamped plate are read from a 
data file and estimated modal frequencies and loss factors for various layer 
thicknesses are output to another file. Material properties of ISD— 112 are 
computed using University of Dayton data and curve-fitting equations to the 
reduced frequency nomogram [Ref. 9,12]. The units used in this program 
are pounds, inches, and seconds. 
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on o onnoonoonooooooooonooonnononnnrinnooooooooooooononnnr 



pRoonAii cARrri 

t xi(M *;()(>!xy?t)(<x><)nix)n>K)>*KX*xx)«K 

C IMIS PROORAn IS TO rAU'ULAIP. 5V51GI1 I OSS PACIORS ( OR VARIOUS 
C I UU’GRA I URCS AMI) IIUDIS UP A COIIS I R A I Mi; I) LAYI.R VISCOGLASIIC 
f. DAMI'IIIG SYSmi MSlMu IMG R US S - K GRl 1 1 M - UMOA R GOUAIIOMS. S'iSIfH LOSS 
C (ACIORS ARG COni'UIPL) I OR A IG DP L R A I MR G RAMOG OP SO 100 DGGRGGS 
'• (ARGMMGII. DASG PLAIC IMtCKMGSSCS VARY PROD O.S75 10 0.75 IMCMPS 

III 0.J25 IMCRGDGMIS. I IM< P.ACII DASE PLAIC TIIICKMGSS, 1IIE VISCOGLASllC 
1II1CKIIGS5 IS VARICO PROD 0.0)5 10 0.U60 IIICIIGS III 0.015 IIICII 
JMCRGDCIIIS. lOR GACII VISCOGLASllC llllCKIIISj IMP; COM 5 I R A 1 II 1 IIG 
LAYGR IIIICKIIESS IS VARIED PROD 0.0625 10 0.25 IIICIIES III 
0.0625 IIICII I IICRGI1GII I S . 

IIIIS PROGRAIt APLIES 10 A rKCE-PREE-rREE-rREE PLATE. 

THE VISCOELASTIC MATERIAL IS 311 ISD-1I2. 

THE FOLLOMIIKI COPITICl PUTS ARE DEPIIIED: 

El - YOllllG'S MODULUS OP DASE PLATE (PSD 
E2 = YOUIIG'S MODMIllS OF VISCOELASTIC LAYER (PSD 

E3 = YOUIIG’S MODULUS OF COIIS T RA I II 1 110 LAYER (PSD 

02 = SHEAR MODULUS OF VISCOELASTIC MATERIAL (PSD 
Mill = POISSOII’S RATIO OF DASE PLAIE 
IIU2 - POISSOII’S RAIIO OF VISCOELASTIC MATERIAL 
HI - IIIICKIIESS OF DASE PLATE (III) 

112 = IIIICKIIESS OP VISCOELASTIC LAYER (III) 

113 = THICKIIESS OF COIIS I K A I II I IIG LAYER (III) 

HIOT = TOTAL PLATE THICKIIESS (111) 

RIIOl = DEIISITY OP BASE PLAIE ( L D P-S ECx x 2/ I II* x A ) 

RI102 = DEH311Y OF VISCOELASTIC MATERIAL 
KII03 = DEIISITY OF COIIS TR A I II I IIG LAYER 

T 0 , FROM, IIROH, II, ML , ETFROL , SL , SH, FROL , 8 C ARE COEFP ICI EHT S 
FOR THE REDUCED FREQUEHCY HOMOGRAM EQUATIOKS 
FP = MODAL FREQUEHCY OF THE UIIDAMPED BASF. PLATE (HZ) 

FCP = MODAL FREQUEHCY OF THE COMPOSIIE PLATE (HZ) 

ETA2 = LOSS FACTOR OF VISCOELASTIC MAIERIAL 
ElAS = SYSTEM LOSS FACTOR 
KQR =.HAVE HULinPR OF PLATE 
GC = GRAVITATIOHAL COHSTAHT (III/SEC**2) 

T = TEMPERATURE OF VISCOELASTIC MATERIAL (DEG F) 

HP = MODAL FREQUEHCY OF UIIDAMPED PLATE (RAD/SEC) 

FP = MODAL FREQUEHCY OF UIIDAMPED PLATE (HZ) 

MCP = MODAL FREQUEHCY OF DAMPED PLATE (RAD/SEC) 

FCP = MODAL FREQUEHCY OF DAMPED PLATE (HZ) 

THE UlllTS USED IH THIS PROGRAM ARE LB, IHCH, SEC, DEGREES F 

UIIDAMPED MODAL FREQUEIICIES ARE FROM A FIIIITE ELEMEIIT AHALYSIS 
OF THE FREE-FREE-FREE-FREE PLATE. 

RESULTS ARE OUTPUT TO DATA FILE "LOSFCTR DATA" 

PRIOR TO RUHIIIIIG THE PROGRAM TYPE THE COMMAHD 
"FILEDEF LOSFCTR DISK LOSFCTR DATA" 

FIIIITE EIEMEIIT MODAL IREQUEHCIES ARE IIIPUT FROM FILE. 

"PLTFRQ DATA" 



K**X**X.*XX**,*XXXXKXXKKX*X*XXKX*K*XKKX*X*****XKXX**XXXX.X**X* 

REAL El , E2, F;3, HUI ,IIU2, HI , 112, 115, RIIOl , RH02, RH03, III OT 

REAL 1 0, FROM,MROM,H,ML, EIFROL,SL , SH, FROL ,C 

REAL FP, FCP, ETA2, El AS, KQR, P I, GC,T 

REAL FRI0,E1A210,SUII1,SU1)2,SUD5,M10 

real M2I , 1151 , G, Cl , AL pure, AI.PIIIH, BRE, Bill, del re 

REAL DELHI, EHCUDE,5UDQ, SUD5,ITP,1ICP, DEHS, SUD6, SUB7,SUB8 

IHIEGER V 

DIMEIISIOII HP(8) 

PI =<i . xaT AIK I . ) 

OPCH(IIH1I = IO,F1IC='LOSFC1R',STATUS=’01D') 
OPEII(UHIl-lI,riLC=’PLirRQ',STAlUS=’OLD’) 

ASSIGH M.MGRIAL COHSIAIIIS FOR BASE P L A T E , COIIS T R A 1 II I IIG LAVER 
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C Alin VISCOELASTIC 
C 

n =1 . 0E7 
Ml ^0 . 7 so 
RM01 = .0068 
IIUI ^ . SS 
RM02= .035 
ES = 1 .0E7 
RII05= .0968 
GC=386 . 

IIU2=.5 

C 

C DCrillE COfISTAIITS FOR MOHOGRAH CQUATIOIIS 
C 

T0 = 1 0^ . 0 
rROM = 2 . 

MR0M=688 . 96 
11^ .275 
ML =8 .7 
EirR0L=l . 08 
SL=0.65 
SH=-.55 
FROL =5000.0 
C = 2.5 
C 

C BASE PLATE TMICKIIESS LOOP 
C 

H1=0 . 375 
no 100 J = 1 , 6 

HRITEC 10,703 ) * BASE PLATE TMICKIIESS, = 

703 rORMAT(A27,F6,3) 

HR I TEC 10, X) 

C 

C VISCOELASTIC THICKNESS LOOP 
C 

M2= 0.015 
DO 150 V=l,6 
C 

C rONSTRAINIllG LAYER THICKNESS LOOP 
C 

H3=0.0625 
DO 200 L=l,6 

IIRITEC 10, 702) ’VISCOELASTIC LAYER THICKNESS, H2=*,H2 
HRITEC 10,?<) 

IIRI lEC 10,702) ’CONSTRAINING LAYER THICKNESS, H3=’,H3 
702 FORMAT(A35,F6.6) 

HRITE(10,>^) 

IIRITEC 1 0, 7 00 ) ’TEMP’, ’MODE’, ’FCP’, ’ETA2’, ’ETAS’, ’G2* 

7 00 FORMAT CAS, 3X, A5,6A15) 

C 

C PLATE MODE LOOP, READ MODAL FREQUENCY AND COMPUTE HAVE NUMBER 
C 

HT0T = HHH2MI3 
DO 300 1=1,8 
READC 1 1 , M) l|p( I ) 

C 

C CALCULATE MODAL FREQUENCY AND HAVE NUMBER OF BASE PLATE 
C . 

SUB7 = SORTC C E1kCHT0Txk3)kGC)/C 1 2. k( 1 , - flUl ) mrmoI mHI ) ) 

KQR = MPC I )/SUB7 
FP=HPC D/C2 ) 

C 

C TEMPERATURE LOOP 
C 

T = 30 . 0 

DO 600 K=1 , 15 
C 

C CALCULATE PROPERTIES OF VISCOELASTIC FOR GIVEN TEMPERATURE AND 
C MODE 
C 

501 FR10 = LOG10CFP)-C12.?<CT-T0))/C 52 5.-i!-T0) 

FR = 10 . .^KU RIO) 
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A- C I R1 0-U)Ol 0( r ROL ) )/C 

SUni-C?‘{(GUjH)?< A KSl-r.li))^(l.“3QRrCl.»AXK2. )))/?. 

Cl A21 U = L C)(n 0 ( Cl TROL ) * SUbl 

c 

C Vi:.C0ELA31 1C loss I'AClUR 

c 

CIA2^10. )txcriA210) 

SUB2 = 2 . JfLOGHUnRUI1/HL) 

SUD5=1 . » Cl R(IM/rR)?<^II 
niO=LOGl OCML )43UD2/Sun3 
C 

C VJ5COELA51IC 5IICAR MOIIULUS 
C 

G2 = lO.K)((niO) 

c 

C VISCOELASTIC YOUIIG'S MODULUS 

n 

E2 = G2k2.x( 1 .HIIU2) 

C 

C CALCULATIONS FOR DAMPED PLATE AT TEMPERATURE T USING ROSS- 
C KCRIIlll-UIIGAR EQUATIONS 
C 

G = G2/( r3xNJ)(H2XKQR) 

U21 = (HHU2)/2. 

I!51=N2+CI11 ni3)/2. 

Cl =C1XN1X( 1 . -»G)4GxE3xN3 
D = GXEIA2X{ F1XHME3XN3) 

AIPIIRE = GXE1MI1XE3X||3X(I131XK2.)><(C1 + DXETA2) 

ALPNIM = GXElXiri XE3XH3XOI31XK2 . )xClxEl A2 
DRE=E1HN1XE2XN2XH31X(C1 +DXETA2) 
BIM=E1XM1XE2XN2XH51X(C1XETA2“D) 
suns = 2 . XGXE2XI12XE3XM3XN21 XN31 XC 
DELRE-SUB8XC C I .-( ETA2XX2 . ) ) + ( DX(2 . XETA2) ) ) 

DrLIM = SUB8x(Clx( 2 . XETA2)-C DX{ 1 . -{ ETA2XX2 . ) ) ) ) 

SUDA = ( 12 ./(Cl XX2. -* DXX2 . ) ) x ( A L PHR E- BRE- DE L RE ) 

EHCUDE=C E1XCM1XK3 . ) ) ^ E3x ( M3xx 3 . ) 45UB<^ 

C 

C MODAL FREQUENCY OF DAMPED PLATE 
C 

DENS=RM01xH14RN02xN2^M3xRN05 

SUBS^C EHCUBExgC)/(12.x( 1 . -NUl xx2 . ) x DENS ) 

HCP=KQRxSQR 1 (SUDS) 
rCP=HCP/(2.xPl ) 

C 

C COMPARISON OF FP AND TCP 
C 

IF (ABS( 1 .-FP/rCP) .LE. 0 . 1 ) THEN 
GOTO 500 
ELSE 

FP=FCP 
GOTO SOI 
ENDIF 
C 

C COMPUTE SYSTEM LOSS FACTOR 
C 

500 SUD6 = (12./(Clxx2. + nxx2. ) ) X ( ALPN I M-D I M- DEL I M ) 

El AS=( 1 ./ENCUBE)XSUD6 
C 

C PRINT RESULIS 
C 

HR 1 TEC 10, 7 01 ) T, I ,FCr, ETA2, ETAS,G2 
701 FORMAT (5X,F7 . 3 , 2X , 1 2 , 3X , N E 1 5 . A ) 

C 

C NEXT lEMPERATURE 
C 

T = T + 5 . n 

^iOO CONTINUE 
C 

C NEXT MODE 
C 

HRITEC 10, X) 

300 CONTINUE 
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RCiiuiDcuni r = i 1 ) 

c 

C IltXT C0U31RA1MIIIG LAYLR llllCKIir.33 
C 

= n . 0625 

200 comiliur: 
c 

C IlEXr VI3C0ELAS11C LAYER TIIICKNC3S 
C 

112 = 112»n.0l5 
150 COlMlIiUC 
C 

C IICXT BASE PLATE 1IUCKIIE3S 
C 

Ml =HI < n . 1 25 
100 COMTIMUE 

Cl OSE(IM{1T = 10) 

CLOSECUni T=1 1 ) 

EMU 
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APPENDIX B 



FORTRAN PROGRAM USED TO COMPUTE MODAL LOSS FACTORS 
FOR THE DOUBLE DAMPING LAYER CONFIGURATION DESIGN 

This program uses the Ross-Kerwin-Ungar equations of Chapter 3 
to compute estimated modal frequencies and loss factors for the double 
layer configuration. Modal frequencies for an undamped plate are read 
from a data file and estimated modal frequencies and loss factors for 
various layer thickness combinations in the double layer configuration are 
output to another data file. Material properties of ISD-112 are computed 
using University of Dayton data and curve-fitting equations to the reduced 
frequency nomogram [Ref. 9,12]. The units used in this program are 
pounds, inches, and seconds. 



PRUGUAM IllOl YR 
C 

c 

C nU3 PROGRAM COMPiJTCS THE SYSTEM I OSS PAnOR AllO MOUAl TREQUrilCY Of 

C A nOUniF. COnSlRAIHEI) LAYER V I SCO EL A S 1 I C AL L Y nAMPEI) PLAIE. lilt LOSS 

C TACMURS APE COMPUICD EUR A SPECIEJC P L A I E/ DAI IP I IIG LAYER 

C conn r.URA 1 I Oil AHD over a TEMPLRAIUR'E range 01 50-100 DEGREES 

C EAHEIIIIEiT. 

C 

C nils PROGRAM APPLIES TO A E R E E - E R E E - E R E E - E R E E I'LAIE AMD IDE UHDAMPED 

C MODAL IREQUEIICIES ARE DE 1 E Rf 11 II E: I) EROII EIIIIIE ELEMCIIl AllALYSIS AIIU 

C ARE READ INTO THIS PROGRAM EROM FILE ’PLlfRQ’. 

C 

C THE VISCOELASTIC MATERIAL IS 3M ISD-112. VISCOELASTIC MATERIAL 
C DAIA IS EROM UIIIVERSIIY OF DAYIOII RESEARCH IHSHIUIE 
C 

C THE UlllTS USED III TtllS PROGRAM ARE LD, IHCfM SEC, AHD DEGREES f. 

C 

C LOSS FACTORS, DAMPED PLATE MODAL FREQUENCIES, AHD ISD-112 PROPERTIES 

C ARE COMPUTED 111 SUDROUTIHE *RKU* FOR EACH ll-TII CONS I RA I II E D LAYER 

C SYSTEM. UPON COMPLETIOII OF EACH TEMPERATURE COMPUTATION THE SYSIEM 
C LOSS FACTOR AHD CORRES PON D I HG DAMPED PLAIE FREQUENCY ARE MRITTEN TO 
C FILE MIIOLYR DATA* . 

C 

C THE rOILOMIllG COErriCIENTS ARE DEEIMED: 

C El = YOUNG’S MODULUS OF ITASE PLAIE (PSI) 

C E5 = EQUIVALENT tOUIIG'S MUDUIUS OF N-TIl CONSTRAINED LAYER SYSTEM 

C E5PRI1 = YOUNG'S MODULUS OF THE CONS T R A I H I IIG LAYER IN ll-TH LAYER 

C ENCODE = EQUIVALENT STIFFNESS OF CONSTRAINED LAYER SYSTEM 

C ETAS = SYSTEM LOSS FACTOR 

C ETAS = LOSS FACTOR OF ll-TN CONSTRAIIIED LAYER SYSTEM 

C ETA5PM = LOSS FACTOR OF CONS T RA I II I NG LAYER IN THE N-TN LAYER ( = 0 

C FCP = FREQUENCY OF THE DAMPED PLATE (HZ) 

C EP = FREQUENCY OF UNDAMPED PLATE (HZ) 

C GC - GRAVITATIONAL CONSTANT 

C 111 = THICKNESS OF BASE PLATE 

C H2 = THICKNESS OF 1ST VISCOELASTIC LAYER 

C 115 = THICKNESS OF N-TH CONSTRAINED LAYER SYSTEM 

C IKi = TOTAL DAMPED PLATE THICKNESS 

C HIPRM = THICKNESS OF BASE LAYER IN N-TM LAYER 

C M2PRM ^ THICKNESS OF VEM IN N-TH LAYER 

C N5PRM = THICKNESS OF CONSTRAINING LAYER IN N-TH LAYER 

C KQR = HAVE NUMBER OF UNDAMPED PI ATE 

C NUl = POISSON'S RATIO OF RASE PLATE AHD CONST RA 1 II 1 IIG LA\CRS 

C NN2 = POISSON'S RATIO OF VISCOELASTIC MATERIAL 

C IIU2PRM = POISSON'S RATIO OF VEM IN ll-TII LAYER 

C RHOl = DENSITY OF BASE PLATE AND CONS 1 HA I II I NG LAYERS 

C RH02 = DENSITY OF VISCOELASTIC MATERIAL 

C RN02PM = DENSITY OF VEM IN N-TH LAYER 

C RN05 = DENSITY OF N-TH LAYER (= RIIOl ) 

C RII05PM = DENSITY OF N-TH LAYER CONSTRAINING LAYER 

C T = TEMPERATURE VARIABLE 

C HP = FREQUENCY OF UNDAMPED PLATE (RAD/SEC) 

C 

C C, ETFROl , FROL , FROM, ML , MROII, II, SH, SI , T 0 = COEFFICIENIS FOR THE 

C REDUCED ERCQUENCY NOMOGRAM EQUATIONS 

C 

c 

REAL C, El , E3, E5PRM, EHCUD E , E T AS , E T A 5 , E T A5PM , E T FROL , ECP, EP, FROL 
REAL FROM, GC, II I , H2 , II 5 , H A , N 1 PRM , II2PRM , II 5f'RM, K OR , ML , DRUM, II, NUl 
REAL IIU2, HU2PM, P I , RllOl , RN02 , RHOZril, Sll , SL , SUD 1 , T , T 0 , 1 IP 
REAL RN05,RHU5PM 
C 

DIMENSION MP( 17 ) 

PI =A . OxATAIK 1.0) 

C 

0PEN(NHIT = 10,FII. E='PLTERQ' ,STATNS='OI D' ) 
0PEN(UN1T=11,F1LE=*IM0LYR',STA[US='0LD') 

C 

C ASSIGN LAYER i II I CKNESSE5 
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c 



111 =0 . 5 
ll?=U 

0 .09^7 5 
I12I’K11 = 0 .015 
li:i'KM=0 . U9575 
ii: = m I'Knui^iMun n jI'Rii 
IKi = ltl < 112 • 113 



c 

C A5SIGII flATCRIAL COIlS I All 1 S 
C 

E) . 0E7 
RllUl =0 . 0968 

mil = 0 . 35 

R1102 = 0 . 035 

11112 = 0. <49 

R1102rVl = 0 .035 
HU2Pil=0 , 5 
RI103 = 0 . 0968 
E3PRM= 1 . 0E7 
Rll03ni=0 . 0968 
E7 A3f'H=0 . 0 
ETA3 = 0.0 
C 

GC = 386 . 0 
C 

C DEriflE VISCOELASTIC COIISTAHTS FOR IIOMOGRAM EQUATIONS 
C 

T0= 106.0 

ruon= 2 . OCA 

MRUn=688 .96 
11 = 0.275 
ia=8 . 7 
ETFROL = 1.08 



SL = 0 . 65 
SII = -0 . 55 
FROL=5000.0 
C = 2 .5 



C 

C HRITE PLATE CHARACTERISTICS 
C 



700 

701 



HRITCC 11,700) 
HR1TE(11,700) 
mu TEU 1 , 7 0 0 ) 
IIRnE(ll, 7 00 ) 

IIRI I E( 1 1 , 7 00 ) 

IIRJTEC 1 1 , X) 
mUTEC 1 1 , X) 
rORllAT(A12,F6 .5) 
11RITE( 11,701) * lEMP ' , 
FORIIAK A8, 3X,A5,2A15) 
IIRITEC 11, X) 



MU =MI1 
M12 =* ,H2 
UllPRlHE =V 
M12PRIIIE = * 
M13PRJIIE =* 



MODE 



HIPRM 

M2PRI1 

HSPRII 



* , *rcp 



C 

C PLATE MODE LOOP 
C 



* El AS* 



DO 100 1 = 1, 17 
READ( 10, X) HP( I ) 

C 

C CALCULATE MODAL FREQUEIICY AND HAVE HUMBER OF UIIDAMI’EU PIAIE 
C 



SUB1=S0RT( (F:1x(I16xx 3))<GC)/( 12.X( 1 . -IIUI xx2 ) xRHO 1 xl|6 ) ) 
KQR=HP( I )/SUDl 
rP = HP( I )/(2 . xpl ) 

C 

C TEMPERATURE LOOP 
C 



T = 30 .0 

DO 200 K= 1 , 1 5 
C 

C COMPUIC SYSIEM LOSS FACTOR AND SIIFFIIESS FOU N-llI LAYER 
C 



CALL RKU( 111 PRIM M2PRi:, M3PRM, LI , RllOi 



iiui , m;2PM, Riio2ni, E3PRM, RiiUiPn, 



1 02 



Cl 1 A.M'M, I , I r , KgK, C 1 AS, LIICUDt: , rci\ 1 0 , IROIU IIKUn, IMIL , t H ROL , 5L , Sll, 

Cf ICUl , C ) 

c 

C COlIVLRI RCSULIS IROIl II-IM LAYER CALCULAIIOII TO TOTAL PIAIE 
C 

n As= n . 0 

1 1* = lire 1 )/( 2 . xpi ) 
c ir=rcr 

rs^i.ncunE/(iisxK3) 

c 

C COMPUTE SYSTEM LOSS TACTOR AIID ERCQUEIICY TOR lOIAL PLATE 
C 

CALI. RKU( Ml , 112, M3, El , Rlini , MUl , MU2, RII0 2, E .’i , Rl 10 3 , E 1 A 3 , T , EP , KQ R , E T A 5 , 
CGMCUDE, TCP, TO, T ROM, IIROM, II , ML , E T F ROL , SL , SI I , f PvOL , C ) 

C 

C PRIIIT RESULTS 
C 

»'RITE(Il,7U2) T, I, rep, ETAS 
702 r0RMAT(5X,F7 .3,2X;I2,3X,2E15.^) 

C 

C NEXT TEMPERATLIRE 
C 

T=Tf5.0 

200 COlinilUE 
C 

C IIEXT MODE 
C 

MRI TEC 1 1 , X ) 

100 COIITIIIUE 
C 

CLOSE(UMIT=I 1 ) 

CL0SE(UriIT = 10) 

EIID 

C 

C 

SUDROUTIME RKU ( H 1 , M2 , II 5 , E 1 , RHO 1 , IIU 1 , NU2 , RII02 , E5 , RttO 3 , ETA3. T,FP, 
CKQR, ETAS, EMCUBE, FCP , T 0 , F ROM, MROM, M, ML , E T EROL , SL , SM , F ROL , C ) 

C 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

c 

C THIS SUBROUFIME CALCULATES VISCOELASTIC PROPERTIES BASED Oil HIE 
C UNIVERSITY OF DAYTON REDUCED FREQUENCY NOMOGRAM EQUATIONS, AND THfN 

C CAtCULATES PLATE STIFFNESSES AND LOSS FACTORS BASED Oil THE ROSS- 

C KERMIN-UNGAR EQUATIONS. 

C 

C THE FOLLOHING ADDITIONAL VARIABLES ARE DEFINED FOR USE IN INIS 
C SUBROUTINE: 

C 

C A = COEFFICIENT FOR NOMOGRAM EQUATIONS 

C ALPIIIM, ALPHRE = IMAGINARY AND REAL COMPONENTS OF COEFFICIENT 

C ALPHA IN THE RKU EQUAriONS 

C DIM,DRE = IMAGINARY AND REAL COMPONENTS OF COEFFICIENT 'D' IN 

C THE RKU EQUATIONS 

C Cl,D = COEFFICIEIMS FOR RKU EQUATIONS 

C DEL1M,DELRE = IMAGINARY AND REAL COMPONENTS OF COEFFICIENT 

C DELTA IN THE RKU EQUATIONS 

C DENS = CminiNATION OF MATERIAL DENSITIES USED TO COMPUTE THE 

C FREQUENCY OF THE DAMPED PLATE AND DAIIPHIG LAYERS 

C ENCUBE = EQUIVALENT STIFFNESS OF DAMPED PLATE AS COMPUTED 

C USING RKU EQUATIONS 

C ETA2 = LOSS FACTOR OF VEM COMPUTED IN NOMOGRAM EQUATIONS 

C ETA210 = LOGIOCEIA?) 

C EIAS = SYSTEM LOSS FACTOR COMPUTED BY RKU EQUATIONS 

C FCP = MODAL FREQUENCY OF DAMPED PLATE ( N7 ) 

C FP = MODAL FREQUENCY OF UNDAMPED PLAIE (N2) 

C FR = REDUCED FREQUENCY OF VEM 

C FRIO = LOGin(FR) 

C G2 = SHEAR MODULUS OF VEM 

C 1121, H3I = RKU EQUAIION fOEFFICIENTS 

C MIO = L0G1U(G2> AS COMPUIED BY NOMOGRAM. EQUATIONS 
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c iic(’ = rRLQurnCY or DAiirci) riAif: (rad/sco 

c 

c MXK)(KXX)()‘KXX?<)(X?<XXXKKX)n<MKXKXXKX?(XXXXX><KXK?<XKKKKKKKXKKXKK)()(K)(KKKKMK)(K 

c 

REAL A, Al rmiu ALTIIRE, DIM, DRE , C , C 1 / D . DLL I M, DEI RE , OEMS, El , E2, E5 
REAL EIICUnC,ETAe,LlA3,r:lA210,ElAS,r:n R()L,rCP,rP,rR, rR10,rRl)L,rROI1 
REAL G. 02. (>c jn . 112, 113. 1121 . M3 1 , KQR . tU 0 . IIL . IIROH, N. MUl . MU2. PI , RMOl 
REAL RII02, RM0 5, SM. SI , 5U fi 1 , SU D 2 . 5U 11 3 , SU D'l , SU D 5 . 5U D6 . SU D 7 , SU D8 
REAL SUB9, sum 0. 5UU1 1 , T.lO.llCP 
C 

PI=A . OKATAMC 1.0) 

GC=386 . 0 
C 

C CALCULATE PROPERTIES OF VEM FOR GIVEN TEMPERATURE AND MODE 
C 

501 FRl 0 = L0G1 0( FP )-( 12. x(T~T0))/( 525. -tT-TO) 

FR=10 , KXCFRl 0) 

A-:(f R10 -l0G10(FR0L ) )/C 

SUm=Cx((SL4SM)xAMSL-3H)N(l.~sgRT(1.4Axx2)))/2. 

ETA210 = LOG10( ETFROL ) + SUDl 
C 

C VISCOELASTIC LOSS FACTOR 
C 

ETA2=10.xx(ETA210) 

C 

C VISCOELASTIC SMEAR MODULUS 
C 

SUD2=2.XLOG10(MROM/ML ) 

SUB3=1 . + ( FROM/FR)XXM 
M10=LOG10(ML )+SUD2/SUD3 
G2=10 .xx(MlO) 

C 

C VISCOELASTIC YOUNG'S MODULUS 
C 

E2=G2x2.x(1.4MU2) 

C 

C CALCULATIONS FOR DAMPED PLATE USING RKU EQUATIONS 
C 

G=G2/(E3XH3XH2XKQR) 

M21 = (HHH2)/2. 

H51=H2 + (HHH3)/2. 

C1=E1XM1X(1 .+G)+GXE3xH3x(1 .-ETA2XETA3) 

D=GxE1xM1xETA2+GXE3xM3x( ETA2+ETA3) 

SUDA=GxE1xM1xE3xH3X(M31 xx2) 

ALPHRE=SUDAX(Cl x( 1 . -ETA2XET A3 )+Dx ( ETA2+ ETA3 ) ) 

ALPHIM = SUB^k(C1x(ETA2^ ETA3) ) 

SUB5=E1xM1KE2xM2xM31 

BRE=SUD5X(C1+DXETA2) 

BIM=SUD5x(C1xETA2-D) 

SUB6=2 . XGXE2XM2XE3XM3XH21XM31 
SUD7=1.-2.xeTA2xETA3-(ETA2XX2) 

SUB8 = 2 .XETA2 + ETA3-C ETA2XX2)XEIA3 
DELRE = SUD6X(Cl XSUB7 + DXSUB8 ) 

DEL IM=SUB6X(C1XSUD8- DXSUD7 ) 

C 

SUD9=(12./(C1xX21Dxx2))X(ALPMRE-BRE-DELRE) 

EMCUDE = E1X(H1XX3)^ E3x(H3xx3)4SUB9 
C 

C MODAL FREQUENCY OF DAMPED PLATE 
C 

DENS=RM01XM1+RM02xM2+RM03xM3 
SUB10=( EMCUDExGC)/( 1 2 , X ( MU 1 X x 2 ) x DENS ) 

MCP=KQRXSQRT(SUD10) 

FCP=HCP/(2.XPI) 

C 

C COMPARISON OF FP AND FCP 
C 

IF(ADS( 1 .-FP/FCP) .LE. 0.10) TMEN 
GOlO 500 

ELSE 

FP=rCP 
GUiO 501 
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non 



UflDIF 



COrirUlE 5Y5TEf1 L0S5 TACIOR 

1300 5UBll=(12./(Clx>(2<U?<x2))x(AlP}nM -DIM -DELI M) 

ETA5 = ( 1 ./EMCUDE )xc (E5><H5><ETA5)»SUD11) 

C 

CflD 
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APPENDIX C 



DESIGN DRAWINGS FOR THE MACHINING 
OF THE FLOATING ELEMENT AND 
POCKET PLATE CONFIGURATIONS 

The drawings shown in Figures C.l and C.2 were used to machine 
the pocket plate and floating element plate used in the experiments. These 
drawings are included to show the relation of the pocket for the ISD— 112 to 
the cover plate and how the viscoelastic was protected from the heat of 
welding. 




Figure C.l. Design Drawing of the Pocket Plate Configuration. 
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1.91 




Figure C.2. Design Drawing of the Floating Element Configuration. 
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APPENDIX D 



REPRESENTATIVE MSC/NASTRAN DATA DECK FOR 
THE DAMPING CONFIGURATIONS 

This data deck was used to compute the modal frequency response of 
the single layer damping configuration and is a representative sample of 
the NASTRAN decks used for the other finite element models. The values 
in the damping table are from a curve fit to the modal loss factors estimated 
from the modal strain energy method. Since the data deck for the normal 
mode and modal strain energy extraction is virtually identical to this deck, 
the Case Control deck commands for the normal mode analysis are 
included, but are commented out. 

The OUTPUT request provides data for an x-y plot of the modal 
frequency response. 

The units used in this deck are pounds, inches, and seconds. 



Id single, ifr 
sol 30 
$ 

^ tiittiitittiitttittttttttttttttttttttttttttttttttttmttttititmtttttt 

$ THIS DEC! IS TO COMPUTE THE MODAL FIEOUEMCT lESPOISE OF A 15x15 

$ INCH ALUMIIUM PLATE NITH A COISTIAIIED VISCOELASTIC DAMPII6 LATEI. 

$ THE MODEL HAS 252 ELEMEITS NITH 64 ELEMBITS II EACH LATEI. THE QUAD 
$ ELEMEITS AIE OFFSET FIOM THE HEX ELEMBITS AS SUGGESTED BY THE 
$ LITEIATUIE. 

$ 

$ THE MATEIIAL PIOPEITIES OF ISD-112 AIE FIOM THE 3M COlPOlATIOI. 

$ 

$ THE MODAL LOSS FACTOIS II THE DAMPIIG TABLE AIE FIOM A CUIVB-FIT 

$ TO TUB SET OF MODAL LOSS FACTOIS COMPUTED FIOM THE MODAL STIAII EIEIGT 
$ METHOD. 

$ 

$ THE UIITS USED II THIS DECI AlE POUIDS.IICHES, AID SECOIDS 

tlie 3000 
cead 

title ^ MODAL FIEOUEICT lESPOISE / 200 HZ / 3M 

■etOod : 1 

spc - 1 

dload : 10 

frequency ^ 10 

sdaiplng ^ 101 

set 111 ^ 95 

svector ^ all 

acceleratlonfplot, phase) ^ 111 
output (xyplot) 
lyprlnt acce / 95(t3) 

$ 

^ tttttttttttttit 

$ THE FOLLOHIHG LUES AIE THE CASE COITIOL DECK CAIDS FOI THE 
$ lOIMAL MODE EITIACTIOI AID STIAII EIEIGT lEQUEST 
$ letbod : 1 
$ spc : 1 

$ set 10 : all 

$ set 11 = 169, thru, 252 

$ ese : 11 

$ 

^ ttttaaatttttttt 

BEGII BULl 

$ TITLE : SIIGLE LATEI NITH QUAD OFFSET 
$ DATA DECI PIODUCED BY PATIAS VEBSIOI 2.0: 2M0V-89 06:31:48 



GIID 


1 


45.0000 15.0000 0.51500 


GIID 


2 


45.0000 15.0000 0.50000 


GIID 


3 


45.0000 12.5000 0.51500 


GIID 


4 


45.0000 12.5000 0.50000 


GIID 


5 


45.0000 10.0000 0.51500 


GIID 


6 


45.0000 10.0000 0.50000 
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GilD 


7 


45.0000 7.50000 0.51500 


GilD 


8 


45.0000 7.50000 0.50000 


GIID 


9 


45.0000 5.00000 0.51500 


GilD 


10 


45.0000 5.00000 0.50000 


GIID 


11 


45.0000 2.50000 0.51500 


GIID 


12 


45.0000 2.50000 0.50000 


GIID 


13 


45.0000 0. 0.51500 


GIID 


14 


45.0000 0. 0.50000 


GIID 


15 


41.7857 15.0000 0.51500 


GIID 


16 


41.7857 15.0000 0.50000 


GIID 


17 


41.7857 12.5000 0.51500 


GIID 


16 


41.7857 12.5000 0.50000 


GIID 


19 


41.7657 10.0000 0.51500 


GIID 


20 


41.7857 10.0000 0.50000 


GIID 


21 


41.7857 7.50000 0.51500 


GIID 


22 


41.7857 7.50000 0.50000 


GIID 


23 


41.7857 5.00000 0.51500 


GIID 


24 


41.7657 5.00000 0.50000 


GIID 


25 


41.7857 2.50000 0.51500 


GIID 


26 


41.7857 2.50000 0.50000 


GIID 


27 


41.7857 0. 0.51500 


GIID 


26 


41.7857 0. 0.50000 


GIID 


29 


38.5714 15.0000 0.51500 


GIID 


30 


36.5714 15.0000 0.50000 


GIID 


31 


38.5714 12.5000 0.51500 


GIID 


32 


36.5714 12.5000 0.50000 


GIID 


33 


36.5714 10.0000 0.51500 


GIID 


34 


36.5714 10.0000 0.50000 


GIID 


35 


38.5714 7.50000 0.51500 


GIID 


36 


38.5714 7.50000 0.50000 


GIID 


37 


38.5714 5.00000 0.51500 


GIID 


36 


36.5714 5.00000 0.50000 


GIID 


39 


38.5714 2.50000 0.51500 


GIID 


40 


38.5714 2.50000 0.50000 


GIID 


41 


38.5714 0. 0.51500 


GIID 


42 


38.5714 0. 0.50000 


GIID 


43 


35.3571 15.0000 0.51500 


GIID 


44 


35.3571 15.0000 0.50000 


GIID 


45 


35.3571 12.5000 0.51500 


GIID 


46 


35.3571 12.5000 0.50000 


GIID 


47 


35.3571 10.0000 0.51500 


GIID 


46 


35.3571 10.0000 0.50000 


GIID 


49 


35.3571 7.50000 0.51500 


GIID 


50 


35.3571 7.50000 0.50000 


GIID 


51 


35.3571 5.00000 0.51500 


GIID 


52 


35.3571 5.00000 0.50000 


GIID 


53 


35.3571 2.50000 0.51500 


GIID 


54 


35.3571 2.50000 0.50000 


GIID 


55 


35.3571 0. 0.51500 


GIID 


56 


35.3571 0. 0.50000 


GIID 


57 


32.1429 15.0000 0.51500 



1 1 1 



GRID 


56 


32.1429 15.0000 0.50000 


GIIO 


59 


32.1429 12.5000 0.51500 


GRID 


60 


32.1429 12.5000 0.50000 


GRID 


61 


32.1429 10.0000 0.51500 


GRID 


62 


32.1429 10.0000 0.50000 


GRID 


63 


32.1429 7.50000 0.51500 


GRID 


64 


32.1429 7.50000 0.50000 


GRID 


65 


32.1429 5.00000 0.51500 


GRID 


66 


32.1429 5.00000 0.50000 


GRID 


67 


32.1429 2.50000 0.51500 


GRID 


68 


32.1429 2.50000 0.50000 


GRID 


69 


32.1429 0. 0.51500 


GRID 


70 


32.1429 0. 0.50000 


GRID 


71 


28.9286 15.0000 0.51500 


GRID 


72 


28.9286 15.0000 0.50000 


GRID 


73 


26.9286 12.5000 0.51500 


GRID 


74 


28.9286 12.5000 0.50000 


GRID 


75 


28.9286 10.0000 0.51500 


GRID 


76 


28.9286 10.0000 0.50000 


GRID 


77 


28.9286 7.50000 0.51500 


GRID 


78 


28.9286 7.50000 0.50000 


GRID 


79 


28.9266 5.00000 0.51500 


GRID 


80 


26.9286 5.00000 0.50000 


GRID 


61 


28.9266 2.50000 0.51500 


GRID 


82 


28.9286 2.50000 0.50000 


GRID 


83 


28.9266 0. 0.51500 


GRID 


84 


28.9266 0. 0.50000 


GRID 


85 


25.7143 15.0000 0.51500 


GRID 


86 


25.7143 15.0000 0.50000 


GRID 


87 


25.7143 12.5000 0.51500 


GRID 


88 


25.7143 12.5000 0.50000 


GRID 


69 


25.7143 10.0000 0.51500 


GRID 


90 


25.7143 10.0000 0.50000 


GRID 


91 


25.7143 7.50000 0.51500 


GRID 


92 


25.7143 7.50000 0.50000 


GRID 


93 


25.7143 5.00000 0.51500 


GRID 


94 


25.7143 5.00000 0.50000 


GRID 


95 


25.7143 2.50000 0.51500 


GRID 


96 


25.7143 2.50000 0.50000 


GRID 


97 


25.714 3 0 . 0.51500 


GRID 


98 


25.7143 0. 0.50000 


GRID 


99 


22.5000 15.0000 0.51500 


GRID 


100 


22.5000 15.0000 0.50000 


GRID 


101 


22.5000 12.5000 0.51500 


GRID 


102 


22.5000 12.5000 0.50000 


GRID 


103 


22.5000 10.0000 0.51500 


GRID 


104 


22.5000 10.0000 0.50000 


GRID 


105 


22.5000 7.50000 0.51500 


GRID 


106 


22.5000 7.50000 0.50000 


GRID 


107 


22.5000 5.00000 0.51500 


GRID 


108 


22.5000 5.00000 0.50000 
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GEID 


109 


22.5000 2.50000 0.51500 


GIIO 


no 


22.5000 2.50000 0.50000 


GfilD 


111 


22.5000 0. 0.51500 


GfilD 


112 


22.5000 0. 0.50000 


GRID 


113 


19.2657 15.0000 0.51500 


GRID 


114 


19.2657 15.0000 0.50000 


GRID 


115 


19.2857 12.5000 0.51500 


GRID 


116 


19.2657 12.5000 0.50000 


GRID 


117 


19.2657 10.0000 0.51500 


GRID 


lie 


19.2657 10.0000 0.50000 


GRID 


119 


19.2657 7.50000 0.51500 


GRID 


120 


19.2857 7.50000 0.50000 


GRID 


121 


19.2857 5.00000 0.51500 


GRID 


122 


19.2857 5.00000 0.50000 


GRID 


123 


19.2857 2.50000 0.51500 


GRID 


124 


19.2857 2.50000 0.50000 


GRID 


125 


19.2657 0. 0.51500 


GRID 


126 


19.2657 0. 0.50000 


GRID 


127 


16.0714 15.0000 0.51500 


GRID 


126 


16.0714 15.0000 0.50000 


GRID 


129 


16.0714 12.5000 0.51500 


GRID 


130 


16.0714 12.5000 0.50000 


GRID 


131 


16.0714 10.0000 0.51500 


GRID 


132 


16.0714 10.0000 0.50000 


GRID 


133 


16.0714 7.50000 0.51500 


GRID 


134 


16.0714 7.50000 0.50000 


GRID 


135 


16.0714 5.00000 0.51500 


GRID 


136 


16.0714 5.00000 0.50000 


GRID 


137 


16.0714 2.50000 0.51500 


GRID 


136 


16.0714 2.50000 0.50000 


GRID 


139 


16.0714 0. 0.51500 


GRID 


140 


16.0714 0. 0.50000 


GRID 


141 


12.6571 15.0000 0.51500 


GRID 


142 


12.6571 15.0000 0.50000 


GRID 


143 


12.8571 12.5000 0.51500 


GRID 


144 


12.6571 12.5000 0.50000 


GRID 


145 


12.6571 10.0000 0.51500 


GRID 


146 


12.6571 10.0000 0.50000 


GRID 


147 


12.8571 7.50000 0.51500 


GRID 


146 


12.6571 7.50000 0.50000 


GRID 


149 


12.6571 5.00000 0.51500 


GRID 


150 


12.6571 5.00000 0.50000 


GRID 


151 


12.6571 2.50000 0.51500 


GRID 


152 


12.6571 2.50000 0.50000 


GRID 


153 


12.8571 0. 0.51500 


GRID 


154 


12.8571 0. 0.50000 


GRID 


155 


9.64266 15.0000 0.51500 


GRID 


156 


9.64266 15.0000 0.50000 


GRID 


157 


6.42657 15.0000 0.51500 


GRID 


156 


6.42657 15.0000 0.50000 


GRID 


159 


3.21429 15.0000 0.51500 
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GEID 


160 


3.21429 


15.0000 


0.50000 


GIID 


161 


0. 


15.0000 


0.51500 


GEID 


162 


0. 


15.0000 


0.50000 


GEID 


163 


9.64286 


12.5000 


0.51500 


GEID 


161 


9.64286 


12.5000 


0.50000 


GEID 


165 


9.64286 


10.0000 


0.51500 


GEID 


166 


9.64286 


10.0000 


0.50000 


GEID 


167 


9.64286 


7.50000 


0.51500 


GEID 


168 


9.64286 


7.50000 


0.50000 


GEID 


169 


9.64286 


5.00000 


0.51500 


GEID 


no 


9.64286 


5.00000 


0.50000 


GEID 


171 


9.64286 


2.50000 


0.51500 


GEID 


172 


9.64286 


2.50000 


0.50000 


GEID 


173 


9.64286 


0. 


0.51500 


GEID 


174 


9.64286 


0. 


0.50000 


GEID 


175 


6.42857 


12.5000 


0.51500 


GEID 


176 


6.42857 


12.5000 


0.50000 


GEID 


177 


3.21429 


12.5000 


0.51500 


GEID 


178 


3.21429 


12.5000 


0.50000 


GEID 


179 


0. 


12.5000 


0.51500 


GEID 


180 


0. 


12.5000 


0.50000 


GEID 


181 


6.42857 


10.0000 


0.51500 


GEID 


182 


6.42857 


10.0000 


0.50000 


GEID 


183 


6.42857 


7.50000 


0.51500 


GEID 


184 


6.42857 


7.50000 


0.50000 


GEID 


185 


6.42857 


5.00000 


0.51500 


GEID 


186 


6.42857 


2.50000 


0.51500 


GEID 


187 


6.42857 


5.00000 


0.50000 


GEID 


188 


6.42857 


2.50000 


0.50000 


GEID 


189 


6.42857 


0. 


0.51500 


GEID 


190 


6.42857 


0. 


0.50000 


GEID 


191 


3.21429 


10.0000 


0.51500 


GEID 


192 


3.21429 


10.0000 


0.50000 


GEID 


193 


0. 


10.0000 


0.51500 


GEID 


194 


0. 


10.0000 


0.50000 


GEID 


195 


3.21429 


0. 


0.51500 


GEID 


196 


3.21429 


0. 


0.50000 


GEID 


197 


0. 


0. 


0.51500 


GEID 


198 


0. 


0. 


0.50000 


GEID 


199 


3.21429 


7.50000 


0.51500 


GEID 


200 


3.21429 


5.00000 


0.51500 


GEID 


201 


3.21429 


2.50000 


0.51500 


GEID 


202 


3.21429 


7.50000 


0.50000 


GEID 


203 


3.21429 


5.00000 


0.50000 


GEID 


204 


3.21429 


2.50000 


0.50000 


GEID 


205 


0. 


7.50000 


0.51500 


GEID 


206 


0. 


5.00000 


0.51500 


GEID 


207 


0. 


2.50000 


0.51500 


GEID 


208 


0. 


7.50000 


0.50000 


GIID 


209 


0. 


2.50000 


0.50000 


GEID 


210 


0. 


5.00000 


0.50000 
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CQUADl 


1 


1 


198 


CQUADl 


2 


1 


196 


COUADl 


3 


1 


190 


CQUADl 


1 


1 


171 


CQUADl 


5 


1 


151 


COUADl 


6 


1 


110 


COUADl 


7 


1 


126 


COUADl 


8 


1 


112 


COUADl 


9 


1 


98 


COUADl 


10 


1 


81 


CQUADl 


11 


1 


70 


COUADl 


12 


1 


56 


COUADl 


13 


1 


12 


COUADl 


11 


1 


28 


COUADl 


15 


1 


209 


COUADl 


16 


1 


201 


COUADl 


17 


1 


188 


COUADl 


18 


1 


172 


COUADl 


19 


1 


152 


COUADl 


20 


1 


138 


COUADl 


21 


1 


121 


COUADl 


22 


1 


110 


COUADl 


23 


1 


96 


CQUADl 


21 


1 


82 


COUADl 


25 


1 


68 


CQUADl 


26 


1 


54 


CQUADl 


27 


1 


10 


CQUADl 


28 


1 


26 


COUADl 


29 


1 


210 


COUADl 


30 


1 


203 


CQUADl 


31 


1 


187 


CQUADl 


32 


1 


170 


CQUADl 


33 


1 


150 


CQUADl 


31 


1 


136 


CQUADl 


35 


1 


122 


COUADl 


36 


1 


108 


CQUADl 


37 


1 


91 


CQUADl 


38 


1 


80 


COUADl 


39 


1 


66 


CQUADl 


10 


1 


52 


COUADl 


11 


1 


38 


CQUADl 


12 


1 


21 


CQUADl 


13 


1 


208 


CQUADl 


11 


1 


202 


COUADl 


15 


1 


181 


CQUADl 


16 


1 


168 


COUADl 


17 


1 


118 


CQUADl 


18 


1 


134 


COUADl 


19 


1 


120 


CQUADl 


50 


1 


106 


COUADl 


51 


1 


92 



204 


209 


-.25 


CD 

OD 


201 


-.25 


172 


188 


-.25 


152 


172 


-.25 


138 


152 


-.25 


121 


138 


-.25 


110 


121 


-.25 


96 


110 


-.25 


82 


96 


-.25 


68 


82 


-.25 


54 


68 


-.25 


10 


51 


-.25 


26 


10 


-.25 


12 


26 


-.25 


203 


210 


-.25 


187 


203 


-.25 


170 


187 


-.25 


150 


170 


-.25 


136 


150 


-.25 


122 


136 


-.25 


108 


122 


-.25 


91 


108 


-.25 


80 


91 


-.25 


66 


80 


-.25 


52 


66 


-.25 


38 


52 


-.25 


21 


38 


-.25 


10 


21 


-.25 


202 


208 


-.25 


181 


202 


-.25 


168 


181 


-.25 


118 


168 


-.25 


131 


118 


-.25 


120 


131 


-.25 


106 


120 


-.25 


92 


106 


-.25 


78 


92 


-.25 


61 


78 


-.25 


50 


61 


-.25 


36 


50 


-.25 


22 


36 


-.25 


8 


22 


-.25 


192 


191 


-.25 


182 


192 


-.25 


166 


182 


-.25 


116 


166 


-.25 


132 


116 


-.25 


118 


132 


-.25 


101 


118 


-.25 


90 


101 


-.25 


76 


90 


-.25 



196 

190 

171 

154 

110 

126 

112 

98 

81 

70 

56 

12 

28 

11 

201 

188 

172 

152 

138 

121 

110 

96 

82 

68 

51 

10 

26 

12 

203 

187 

170 

150 

136 

122 

108 

94 

80 

66 

52 

38 

21 

10 

202 

181 

168 

118 

131 

120 

106 

92 

78 



CQUAD4 


52 


1 


78 


C0UAD4 


53 


1 


64 


CQUAD4 


54 


1 


50 


CQ0AD4 


55 


1 


36 


C0UAD4 


56 


1 


22 


C0UAD4 


57 


1 


194 


C0UAD4 


58 


1 


192 


C0UAD4 


59 


1 


182 


C0UAD4 


60 


1 


166 


COUAD4 


61 


1 


146 


C0UAD4 


62 


1 


132 


C0UAD4 


63 


1 


118 


C0UAD4 


64 


1 


104 


CQUAD4 


65 


1 


90 


C0UAD4 


66 


1 


76 


CQDAD4 


67 


1 


62 


CQUAD4 


68 


1 


48 


C0UAD4 


69 


1 


34 


C0UAD4 


70 


1 


20 


C0UAD4 


71 


1 


180 


C0UAD4 


72 


1 


178 


C0UAD4 


73 


1 


176 


C0UAD4 


74 


1 


164 


C0UAD4 
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